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Assuming the universe is spatially homogeneous on the largest scales lays the foundation for almost 
all cosmology. This idea is based on the Copernican principle, that we are not at a particularly 
special place in the universe. Surprisingly, this philosophical assumption has yet to be rigorously 
demonstrated independently of the standard paxadigm. This issue has been brought to light by 
cosmological models which cam potentially explain apparent acceleration by spatial inhomogeneity 
rather than dark energy. These models replace the temporal fine tuning associated with A with a 
spatial fine tuning, and so violate the Copernican assumption. While is seems unlikely that such 
models can really give a realistic solution to the dark energy problem, they do reveal how poorly 
constrained radial inhomogeneity actually is. So the bigger issue remains: How do we robustly test 
the Copernican principle independently of dark energy or theory of gravity? 
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I. INTRODUCTION 

The standard model of cosmology is fabulous in its simplicity: based on linear perturbations about a spa- 
tially homogeneous and isotropic background model, it can easily account for just about all observations 
which probe a vast range of scales in space and time with a small handful of parameters. The bigger pic- 
ture which emerges is of a model with an exponential expansion rate for much of the evolution of the uni- 
verse, caused by the inflaton at the beginning and dark energy at the end. We are anthropically selected 
to live in the small era between these phases where structure forms and interesting things happen. Yet 
the physical matter which drives these accelerating periods is not understood at all in a fundamental sense. 
Until they are, the standard model is un- 
fortunately phenomenological in this criti- 
cal aspect. Because of this, the anthropic 
fine tuning seems perverse: the cosmolog- 
ical constant, at odds with its 'expected' 
value by some 120 orders of magnitude, has 
an energy density pA today about the same 
as that of matter p^, despite the fact that 
the ratio of these grows with the volume 
of space: p^l Pm ~ a"^ ~ 1. We are living 
through a phase transition (Fig. 1). Why? 

The problem of understanding the phys- 
ical origin and value of the cosmological 
constant is leading us to reconsider some 
of the foundational aspects of cosmological 
model building more carefully. In particu- 
lar, it is an important fact that, at the mo- 
ment, the spatial homogeneity of the uni- 
verse when smoothed on equality scales ex- 
ists by assumption, and is not yet an ob- 
servationally proven fact established out- 
side the paradigm of the standard model - 
which includes dark energy. Given this un- 
certainty, so-called void models can explain 
the observed distance modulus utilising a 
spatially varying energy density, Hubble rate and curvature on Gpc scales, without any unusual physical fields at late 
times [2-113]. The indication is that models which are homogeneous at early times are incompatible with observa- 
tions, as are adiabatic models with the simplest type of inhomogcneity present at early times [107, 108]. Isocurvature 
degrees of freedom and freedom in the initial power spectrum have not been explored in detail, however, and remain 
possible routes to constructing viable alternatives to ACDM [81, 83, 94]. They are therefore a very significant de- 
parture from the standard model if they can be made to work, and would require a dramatic reworking of standard 
inflation (though there are inflationary models which can produce large spherical structures sec e.g., [114, 115]). 
Irrespective of all this, thej^ have actually been neglected as a serious contender for dark energy because of the anti- 
Copcrnican fine tuning that exists: we have to be within tens of Mpc of the centre of spherical symmetry of the 
backgroimd [21, 25, 61, 84], which implies a spatial coincidence of, roughly, (40 Mpc/15 Gpc)'^ ~ 10~^. This is just 
plain weird. However, it is not hard to imagine a bigger picture where there exist structures as large or larger 
than - our Hubble scale, one of which we are just glimpsing a part of [116]. Perhaps there could be selection effects 
favouring stable solar systems in regions of lower dark matter density (or something), which would normalise the 
spatial coincidence. Who knows? 

While it is still not clear whether these models can really be made a viable alternative to dark energy, these models 
have brought into focus important questions: Is the universe spatially homogeneous and isotropic when viewed on the 
largest scales? Perhaps we really have shown this already? If so, to what level of confidence? 

The applicability of the Friedmann-Lemaitre- Robertson- Walker (FLRW) metric is the underlying axiom from which 
we infer dark energy exists whether in the form of A, exotic matter or as a modification of GR. It is ne(;essary, 
therefore, to try to demonstrate that the FLRW paradigm is correct from a purely observational point of view - without 
assuming it a priori, and preferably independently of the field equations. There are different issues to consider: 

The Copernican Principle: We are not at a special location in the universe. 

The Cosmological Principle: Smoothed on large enough scales the universe is spatially homogeneous and isotropic. 
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FIG. 1: The coincidence problem: why is A as large as it can be? Any larger 
and the de Sitter phase would start before structure forms. (Prom [1].) 
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A textbook formulation of the FLRW metric from these principles starts from the first then uses the high isotropy 
of the Cosmic Microwave Background (CMB), and approximate isotropy of local observables to conclude the second. 
This is correct in an exact sense, as we discuss below: e.g., if all observers measure the distance-redshift relation 
to be exactly isotropic, then the spacetime is exactly FLRW. The realistic case is much more subtle. A statement 
such as 'if most observers find their observables consistent with a small level of anisotropy, then the metric of the 
universe is roughly FLRW when smoothed over a suitably large scale' requires quite a few assumptions about spatial 
gradients which may or may not be realistic, and has only been theoretically argued in the case of the CMB. An 
additional problem here is what we mean by smoothing in a spacetime: smoothing observables is not the same as 
spatial smoothing, and smoothing a geometry is ill-defined and is not the same geometry one arrives at from a 
smoothed energy-momentum tensor (see [117, 118] for recent reviews). We shall not consider this important problem 
further here as it is beyond the scope of this work, and assume that the conventional definition in terms of a smooth 
spacetime is sufficient. 

These important subtleties aside, going from the 
Copernican principle (CP) to the FLRW geometry 
seems reasonable in a wide variety of circumstances, 
as we discuss in detail. Consequently, establishing spa- 
tial homogeneity of the universe - and with it the ex- 
istence of dark energy - can be answered if we can 
observationally test the Copernican assumption. How- 
ever obvious it seems, ^ it is nevertheless a philosophical 
assumption at the heart of cosmology which should be 
demonstrated scientifically where possible [120]. 

The Copernican principle is hard to test on large 
(Gpc) scales simply because we view the imiverse ef- 
fectively from one spacetime event (Fig. 2), although 
it can be tested locally [121]. Compounding this is the 
fact that it's hard to disentangle temporal evolution 
from spatial variation - especially so if we do not have a 
separately testable model for the matter present (dark 
energy!). A nice way to illustrate the difficulty is to 
consider an alternative way of making a large-scale cos- 
mological model. Instead of postulating a model at an 
early time (i.e., an FLRW model with perturbations), 
evolving it forwards and comparing it with observa- 
tions, we can take observations directly as 'initial data' 
on our past lightcone and integrate into the interior to 
reconstruct our past history [45, 54, 80, 91, 122, 123]. 
Would this necessarily yield an FLRW model? What 
observables do we need? Under what assumptions of dark energy and theory of gravity: is a model based on general 
relativity which is free of dark energy a possible solution? While such a scheme is impractical in the near future, 
it is conceptually useful to consider cosmology as an inverse problem, and should be well-posed at least while local 
structure is in the linear regime. 

With these ideas in mind, practical tests of the Copernican assumption can be developed. The are several basic 
ideas. One is to try to directly observe the universe as others see it. An alien civilisation at 2; ~ 1 who had the 
foresight to send us a data file with their cosmological observations would be nice, but failing that placing limits on 
anisotropy around distant clusters can achieve the same ends with less slime. Another is to combine observables such 
that we can see if the data on our past lightcone would conflict with an FLRW model in the interior. This helps 
formulate the Copernican principle as a null hypothesis which can in principle be refuted. A third is to see if the 
thermal history is the same in widely separated regions which can be used to probe homogeneity at early times [124]. 

This review is organised as follows. First we consider what isotropic observations tell us, and review models which 
violate the Copernican principle. Then we discuss general results which help us go from exact isotropy of observables 
to exact homogeneity of space. Finally we summarise the consistency tests available to test the FLRW assumption. 




^ The wisdom of the crowd can give the ctccurate weight of a cow [119]; it would be engaging to see what this method would give us here. 
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II. MODELS WITHOUT HOMOGENEITY AS AN ALTERNATIVE TO DARK ENERGY 
A. Prom isotropic observables to isotropy of space 

Without assuming the Copernican Principle, we have the following isotropy result [122, 125, 126]: 



1 CP 



> spatial isotropy: If one observer comoving with the matter mea- 



our worldline 



[ENMSW] Matter lightcone-isotropy 

sures isotropic area distances, number counts, bulk peculiar velocities, and lensing, in an expanding dust Universe 
with A, then the spacetime is isotropic about the observer's worldline. 

This is an impressive selection of ob- 
servables required. Note that isotropy of 
(bulk) peculiar velocities seen by the ob- 
server is equivalent to vanishing proper 
motions (tranverse velocities) on the ob- 
server's sky. Isotropy of lensing means that 
there is no distortion of images, only mag- 
nification. 

The proof of this result requires a non- 
perturbative approach ~ there is no back- 
ground to perturb around. Since the data 
is given on the past lightcone of the ob- 
server, we need the fully general metric, 
adapted to the past lightcones of the ob- 
server worldline C. We define observa- 
tional coordinates a;° = {w,y,9,(l)), where 
x"" — {6, (f>) are the celestial coordinates, 
w — const are the past light cones on C 
{y = 0), normalised so that w measures 
proper time along C, and y measures dis- 
tance down the light rays {w, 9, (f) = const. 
A convenient choice for y\sy — z (redshift) 
on the lightcone of here-and-now, w = wq, 
and then keep y comoving with matter off 
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FIG. 3: The Copernican principle is hard to test because we are fixed to one 
event in spacetime. We make observations on our past nullcone which sHces 
through spatial surfaces. 



the initial lightcone, so that — 0. (This is rather idealised of course, as redshift may not be monotonic, and caustics 
will form, and so on.) Then the matter 4- velocity and the photon wave- vector are 



U'' = (1 + Z)(1,0,I/'^), ka^daW, l + Z^Uak\ 

where V"" — dx"" /dw are the transverse velocity components on the observer's sky. The metric is 

ds^ = -A^dw'^ + 2Bdwdy + 2Cadx''dw + D^{dn^ + Labdx'^dx'') 

= (l + ^)-2 + 2aF"+5a6^"l/^ B='^, 

dy 



(1) 

(2) 
(3) 



where the expression for A^ follows from Uau" 
of images via the shear of lightrays, 



T; I? is the area distance, and Lab determines the lensing distortion 



Cab 



2B dy ■ 

The number of galaxies in a solid angle df2 and a null distance increment dy is 

diV = Sn{l + z)D^Bdndy , 



(4) 



(5) 



where S is the selection function and n is the number density. 

Before specializing to isotropic observations, we identify how the observations in general and in principle determine 
the geometry of the past light cone w — wq of here-and-now, where y — z: 



Area distances directly determine the metric function D^wq, z, x") 
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• The number counts (given a knowledge of S) determine Bn and thus, assuming a knowledge of the bias, they 
determine B(wo, z, x°')pm{wQ, z, x"), where pm = Pb + Pc is the total matter density. 

• Transverse (proper) motions in principle directly determine V°'{wo,z,x''). 

• Image distortion determines Lab{wo,z,x'^). (The differential lensing matrix aab is determined by Lab,D,B.) 
Then [125, 126]: 

[ENMSW] Lightcone observations =4* spacetime metric: Observations {D, N,V°' , Lab) on the past lightcone 
w = wq determine in principle (gabjU"^, Bpm) on the lightcone. This is exactly the information needed for 
Einstein's equations to determine B, Ca on w = wq, so that the metric and matter are fully determined on the 
lightcone. Finally, the past Cauchy development of this data determines 

9ab ? ) Pm the interior of the past 

lightcone. 

If we assume that observations are isotropic, then 

Momentum conservation and the yy field equation then give the following equations on w = Wq [122, 125]: 

Ca = (I + Z)-^ r{l + z)B,adz (7) 

Jo 



B = ^=2D' 
d^; 



2- / {l + zfDn^dz 
Jo 



(8) 



where a prime denotes d/dz. These imply that B a = = Ca, so that Pm,a = - and hence the metric and matter 
are isotropic on w = wq. This can only happen if the interior of ui = wq is isotropic. If observations remain isotropic 
along C, then the spacetime is isotropic. 



B. Cosmology with spherical symmetry 

Isotropic observations imply spherical symmetry in the presence of dust matter, leading to the Lemaitre-Tolman- 
Bondi (LTB) model, or the 'LLTB' model if we include A. An interesting explanation for the dark energy problem 
in cosmology is one where the underlying geometry of the universe is significantly inhomogeneous on Hubble scales. 
Spacetimes used in this context are usually LTB models - so-called 'void models', first introduced in [2]. These models 
can look like dark energy because we have direct access only to data on our lightcone and so we cannot disentangle 
temporal evolution in the scale factor from radial variations. The main focus has been on aiming to see if these models 
can fit the data without A, thus circumventing the dark energy problem. However, they can equally be used with 
A to place constraints on radial inhomogeneity, though very little work has been done on this [87]. We shall briefly 
review the LTB dust models, as they illustrate the kind of observations required to constrain homogeneity. 

An inhomogeneous void may be modelled as a spherically symmetric dust LTB model with metric 

ds^ = -di^ + -^iM^dr^ + ai {t, rydn^ , (9) 
1 — K(r)H 

where the radial (ay) and angular (aj_) scale factors are related by ay = d{a±r)/dr. The curvature k = K{r) is not 
constant but is instead a free function. The FLRW limit is /t — >■ const., and a± = ay . The two scale factors define two 
Hubble rates: 

H^=H^{t,r) = ^, =if||(i,r)= ^, (10) 

Ctj_ Cly 

The analogue of the Friedmann equation in this space-time is then given by 

- 4 " 4 ' ^^^^ 
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where M = M{r) is another free function of r, and the locally measured energy density is 



8nGp{t, r) 



which obeys the conservation equation 

The acceleration equations in the transverse and radial directions are 



a_L 



M 
'24 



and 



ail 



-Airp - 



M 



We introduce dimensionlcss density parameters for the CDM and 
curvature, by analogy with the FLRW models: 



Hi ■ 



M 

W 



(15) 



using which, the Friedmann equation takes on its familiar form: 
HI 



Hi 

-Lo 



(16) 



so flm{r) + ^^{1") = 1- Integrating the Friedmann equation from 
the time of the big bang ts = tBir) to some later time t yields 
the age of the universe at a given {t, r): 



T{t,r) =t-tB 



a.±{t,r) 



dx 



(17) 



(12) 



(13) 



(14) 




FIG. 4: A void model produced by a Newtonian N- 
body simulation. (From [90].) 



We now have two free functional degrees of freedom: flm{r) and 
tsir), which can be specified as we wish (if the bang time func- 
tion is not constant this represents a decaying mode if one tries 
to approximate the solution as perturbed FLRW [127]). A co- 
ordinate choice which fixes a±{to,r) = 1 then fixes H±g{r) from 
Eq. (17). A value for Hq = H±g{r = 0) is used at this point. 

The LTB model is actually also a Newtonian solution - that is, Newtonian spherically symmetric dust solutions 
have the same equations of motion [128]. This was demonstrated explicitly in an N-body simulation of a void [90] 
(see Fig. 4). 



C. Background observables 

In LTB models, there are several approaches to finding observables such as distances as a function of redshift [5]. 
We refer to [32] for details of the approach we use here. On the past light cone a central observer may write the t, r 
coordinates as functions of z. These functions are determined by the system of differential equations 



dt 1 dr Vl — nr'^ 

dz ^ ^ (l + z)Hii ' dz ~ (l + ^)a|[iJ||' ^ ' 

where if||(t, r) = H^\{t{z),r{z)) = -ff||(-z), etc. The area distance is given by 

dA{z) = a^{t{z),r{z))r{z) (19) 

and the luminosity distance is, as usual dL{z) = (1 + z)^dA{z). The volume element as a function of redshift is given 
by 

dV ^ AndAjz)^ 

dz {l + z)Hn{z)- ^ ' 
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This then imphes the number counts as a function of redshift, provided the bias and mass functions are known; if not, 
there is an important degeneracy between source evolution when trying to use number counts to measure the volume 
element [7, 10]. 

With one free function we can design models that give any distance modulus we like (see e.g., [3, 5, 7, 9, 16, 17, 
19, 20, 22, 32, 39, 41, 43, 60]). In Fig. 5 we show a selection of different models which have been considered recently. 
Generically, those give rise to 'void models': with the bang time function set to zero and we choose r2m(r) to reproduce 




FIG. 5: LTB models have no problem fitting distance data. Left is an attempt to fit the early SNIa data of [129], using a very 
small void with an over-dense shell around it embedded in an EdS model, from [36]. The hope was that we could be located 
in the sort of voids we observe all over the place, giving a jump in the distance modulus which can then fit the SNIa. The 
gap in the data at intermediate redshift was filled by the SDSS SNIa [59] which ruled these out, leaving the possibility of giant 
voids several Gpc across with a Gaussian density profile as an alternative to ACDM (centre, top), with approximate dimensions 
shown (below). If the bang time function is non-zero then the data do not constrain the density to be a void profile (right); 
[55] show that a central over-density can fit the SNIa data of [130]. 



exactly a ACDM D{z), then the LTB model is a void with steep radial density profile, but a Gaussian profile fits the 
SNIa data just as well [60]. 

D. The small-scale CMB & Ho 

The physics of decoupling and line-of-sight effects contribute differently to the CMB, and have different dependency 
on the cosmological model. In sophisticated inhomogeneous models both pre- and post-decoupling effects will play a 
role, but Hubble-scale void models allow an important simplification for calculating the moderate to high £ part of 
the CMB. 

The comoving scale of the voids which closely mimic the ACDM distance modulus are typically O(Gpc). The 
physical size of the sound horizon, which sets the largest scale seen in the pre-decoupling part of the power spectrum, 
is around 150 Mpc redshifted to today. This implies that in any causally connected patch of the Universe prior to 
decoupling, the density gradient is very small. Furthermore, the comoving radius of decoupling is larger than 10 Gpc, 
on which scale the gradient of the void profile is small in the simplest models (or can be by assumption). For example, 
at decoupling the total fractional difference in energy density between the centre of the void and the asymptotic region 
is around 10% [81]; hence, across a causal patch we expect a maximum 1% change in the energy density in the radial 
direction, and much less at the radius of the CMB that we observe for a Gaussian profile. This suggests that before 
decoupling on small scales we can model the universe in disconnected FLRW shells at different radii, with the shell 
of interest located at the distance where we see the CMB. This may be calculated using standard FLRW codes, but 
with the line-of-sight parts corrected for [48, 50]. The calculation for the high-^ spectrum was first presented in [48], 
and further developed in [50, 75, 76, 81]. 
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FIG. 6: Left: The area distance to 2; ~ 1090 in a Gaussian-profiled LTB void model with zero bang time. Adding bumps to 
the density profile changes this figure considerably. Whether the model can fit the CMB lies in the freedom of the value of Hq 
for a measured To and z*. The simplest models require h ~ 0.5. 

Right: The normalised CMB angular power spectrum. The power spectrum is shown against a default flat concordance model 
with zero tilt. There is nothing between the two models for high I, with the maximum difference around 1%. (From [81]) 



For line-of-sight effects, we need to use the full void model. These come in two forms. The simplest effect is via the 
background dynamics, which affects the area distance to the CMB, somewhat similar to a simple dark energy model. 
This is the important effect for the small-scale CMB. The more complicated effect is on the largest scales through the 
Integrated Sachs- Wolfe effect (see [51] for the general formulas in LTB). This requires the solution of the perturbation 
equations presented below, and has not been addressed. 

The CMB parameters (an asterisk denotes decoupling) 



(21) 



are sufficient to characterise the key features of the first three peaks of the CMB [131, 132] (see also [133]). Within 
a standard thermal history, the physics that determines the first three peaks also fixes the details of the damping 
tail [134]. With the exception of dA{z^,), all quantities are local to the surface of the CMB that we observe. 

The parameters given by Eqs. (21) separate the local physics of the CMB from the line-of-sight area distance. These 
parameters can be inverted to provide nearly uncorrelated constraints on (iyi(z*), /b and 77. Specifying asymptotic 
void model parameters to give the measured value of i?* and la/Icq, leaves just the area distance of the CMB to be 
adjusted to fit the CMB shift parameters. This constrains a combination of the void profile and the curvature and 
Hubble rate at the centre. 

A final constraint arises when we integrate out along the past lightcone from the centre out to . In terms of time 
it says that the local time at that z* must equal the time obtained by integrating up along the timelike worldline from 
the big bang up to decoupling. That is, 



dz 
l + z 



nullconc 



(22) 



where x(t, r) = —dt/dr evaluated on the past nuUcone, and i*(z*) is the local time of decoupling at the redshift 
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observed from the centre, which must be equal to 

f°° dT 1 

where the Hubble rate as a function of temperature, H{T), is given locally at early times by 

which also only has dependence on the local parameters r] and fb and with no reference to late times. We have defined 
the dimension-full constants 

n-yh"^ /0.02587V O^/i^ Ofc/i^ 30C(3) 

'^1 = ~ J ' '^v = « 0.227AfeffK7^, vjh = = ^4 m^w^ . (25) 

Note that these have no dependence on any parameters of the model. That is, we are not free to specify the tu's, 
apart from A^eff- These are derived assuming that /b and t] are constant in time. An example from [81] of how closely 
a void model can reproduce the CMB power spectrum found in a concordance model is shown in Fig. 6. 

In [48, 50, 79], it was shown that the CMB can be very restrictive for adiabatic void models (i.e., those with 
J] = const, spatially at early times) when the bang time is zero, the power spectrum has no features, and the universe 
is assumed to evolve from a homogeneous model. We can see this as follows. For a Gaussian profiled void with 
~ 0.1 at the centre, the area distance to the CMB favours a low fi^ asymptotically, or else a low ifg at the centre 
(see Fig. 6). Thus, an asymptotically flat model needs Hq ^ 50kms^^Mpc~^ to get the area distance right. Then, 
if the constraint, Eq. (22), is evaluated either in an LTB model, or by matching on to an FLRW model, it is found 
that the asymptotic value of the density parameter must be high. Thus, in this approximation, we sec that the CMB 
favours models with a very low at the centre to place the CMB far enough away. The difficulty fitting the CMB 
may therefore be considered one of fitting the local value Hq [79], which is quite high. However, [50] showed that 
with a varying bang time, the data for H^^ SNIa and CMB can be simultaneously accommodated. This is because 
the constraint. Eq. (22), must be modified by adding a factor of the diflFcrence between the bang time at the centre 
with the bang time asymptotically, so releasing the key constraint on Hq. 

It was argued in [81] that Eq. (22) can be accommodated by an 0(1) inhomogeneity in the radiation profile, 
T] = ri{r), at decoupling which varies over a similar scale as the matter inhomogeneity - giving an 'isocurvature' void. 
The reason is because the constraint is sensitive to t«/fo ~ 10^^, which is mirrored by a sensitivity in z^, at around 
the 10% level when ^ 1090. Thus [81] argue that a full two-fiuid model is required in order to decisively evaluate 
Eq. (22) in this case and so provide accurate constraints on isocurvature voids, though it remains unknown if this 
gives enough freedom to raise Hq sufficiently. 

An important alternative solution, presented in [46, 94], is to add a bump to the primordial power spectnmi around 
the equality scale. This then allows perfectly acceptable fits to the key set of observables ilo+SNIa+CMB. This is 
particularly important because to produce a void model in the first place, at least one extra scale must be present 
during inflation (or whatever model is used) so it is unreasonable to assume a featureless power spectrum which is 
the case for all other studies. 

Although the simplest models do not appear viable because their local Hq is far too low [48, 50, 79], it is still an 
open question exactly what constraints the small-scale CMB places on a generic void solution. 



E. Scattering of the CMB 

The idea of using the CMB to probe radial inhomogeneity on large scales was initiated in a seminal paper by 
Goodman [135]. In essence, if we have some kind of cosmic mirror with which to view the CMB around distant 
observers we can measure its temperature there and constrain anisotropy of the CMB about them, and so constrain 
the type of radial inhomogeneity prevalent in void models. Observers in a large void typically have a large peculiar 
velocity with respect to the CMB frame, and so see a large dipole moment, as well as higher multipoles; that is, 
observers outside the void will see a large patch of their CMB sky distorted by the void. 

There are several mechanisms by which the temperature of the CMB can be measured at points away from us. The 
key probe is using the Sunyaev-Zel'dovich effect [136, 137], where the CMB photons are inverse Compton scattered 
by hot electrons, typically in clusters. There are two effects from this: it heats up the CMB photons by increasing 
their frequency which distort the CMB spectrum; and it scatters photons into our line of sight which would not 
otherwise have been there. This changes the overall spectrum of the CMB in the direction of the scatterer because the 
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FIG. 7: Off-centre observers typically see an anisotropic CMB sky (top). This causes a large distortion in the CMB spectrum 
(bottom left). The kSZ effect allows us to infer the radial velocities of clusters which have a systematic drift in void models. 
In models with a homogeneous bang time, this is large (middle panel, right, dashed line) and already ruled out by present 
constraints [42, 58]; a small inhomogeneous bang time can drastically reduce this effect (middle panel, right, solid line) without 
changing the late time void (bottom right), and an isocurvature mode can do the same sort of thing [83]. (Figures from; 
left [34]; top right [42]; bottom right [107].) 



primary CMB photons become mixed up with the photons from all over the scatterers' sky. If the scatterer measures 
an anisotropic CMB - so that their CMB temperature is different along different lines of sight - then this distorts 
the initially blackbody spectrum which is scattered into our line of sight. (The sum of two blackbodies at different 
temperature is not a blackbody.) 

The altered spectrum from scattering of CMB photons into our line of sight has two main contributions. The 
contribution from all the anisotropics at a cluster produces a so-called y-distortion to the spectrum [34, 135, 138]. 
Adding up all sources, assuming single scattering results in a distortion proportional to 

yoij^dz^ J d^n'{l + n-n')[^{n,n,z)- ^{n',n,z)f (26) 

where ^(n',n, z) is the CMB temperature anisotropy in the direction n' at the cluster located at z in direction 
n according to the central observer. t{z) is the optical depth. (Note that the thermal SZ effect also produces a y- 
distortion too, but through the monopole temperature rather than temperature anisotropics.) The other contribution 
is from the kinetic SZ effect. This is an effect caused by the bulk motion of a cluster relative to the CMB frame [42, 
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58, 83, 86, 102, 107]. An observer looking in direction n then sees a temperature fluctuation 




Jo 




(27) 



where is the density contrast of electrons along the line of sight. In a void model clusters will have a systematic 
radial velocity Vriz) which can place constraints on the model. In addition, given a primordial power spectrum and a 
way to evolve perturbations, the angular power spectrum associated with a continuum Se{n,z) may be evaluated as 
a correction to the usual CMB Ce's. 

Constraints are impressive. All studies which consider a homogeneous early universe find them severely constrained 
using these effects [34, 42, 58, 82, 86, 102, 135], and effectively rule out such models. However, it was shown in 
[107] that a bang time function with amplitude of order the decoupling time can be used to tune out the kSZ 
signal by fine tuning the peculiar velocity on the past lightcone. Removing the off-centre dipole in this way will 
weaken y-distortion constraints as well. They found, however, that models which did this could not fit the CMB+_ffo 
constraints which requires a Gyr-amplitude bang time in their analysis. (Such a large amplitude bang time induces 
a huge j/-distortion which is ruled out in FLRW [108].) It was argued in [81, 83] that because the temperature of 
decoupling is not in general spatially constant that this should also be used to investigate these constraints, and will 
weaken them considerably; in this interpretation these constraints are really measurements of fb{r) and r]{r). Finally, 
the integrated kSZ constraints [86, 102] rely on structure formation and an unknown radial power spectrum, and so 
have these additional degrees of freedom and problems to consider. 

Nevertheless, scattering of the CMB provides stringent constraints on the simplest voids, and show that for a 
void model to be a viable model of dark energy it will have to have inhomogcneity present at early times as well. 
Furthermore, if one is to tune a model to have vanishing dipole on the past lightcone of the central observer, this will 
typically be possible only for one lightcone, and hence for one instant for the central observer. This will add further 
complications to the fine tuning problems of the models. 



Big-Bang nucleosynthesis (BBN) is the most robust probe of the first instants of the post-inflationary Universe. 
After three minutes, the lightest nuclei (mainly D, ^He, "^He, and ^Li) were synthesised in observationally signiflcant 
abundances [141, 143]. Observations of these abundances provide powerful constraints on the primordial baryon- 
to-photon ratio r] = ni,/n-y, which is constant in time during adiabatic expansion. In the ACDM model, the CMB 
constrains rjcMB = 6.226 ± 0.17 x 10~^° [143] at a redshift z ~ 1100. Observations of high redshift low metallicity 
quasar absorbers tells us D/H= (2.8 ± 0.2) x 10~^ [142] at ~ 3, which in standard BBN leads to rjo = (5.8 ± 
0.3) X 10"^^°, in good agreement with the CMB constraint. In contrast to these distant measurements at z ^ 10^ 
and z ^ 3, primordial abundances at z = arc either very uncertain (D and ^He), not a very sensitive baryometer 
(^He), or, most importantly, in signiflcant disagreement with these measurements - ^Li. To probe the BBN yield of 
''Li, observations have concentrated on old metal-poor stars in the Galactic halo or in Galactic globular clusters. The 
ratio between rj^i derived from ''Li at z = and rjo derived from £) at 2; ~ 3 is found to be "qD/rjn ~ 1.5. Within 
the standard model of cosmology, this anomalously low value for iju disagrees with the CMB derived value by up to 



A local value of ~ 4 — 5 x 10"^*^ is consistent with all the measurements of primordial abundances at z = 0, 
however (see top left panel in Fig. 8). The disagreement with high- redshift CMB and D data (probing rj at large 
distances) shows up only when rj is assumed to be homogeneous on super-Hubble scales at BBN, as in standard 
cosmology. An inhomogeneous radial proflle for r] can thus solve the ''Li problem, shown in Fig. 8 [75] . 



During the process of recombination, when Compton scattering of the electrons and photons becomes low enough 
the baryons become freed from the photons. This 'drag epoch', when the baryons are no longer dragged by the 
photons, happens at a temperature T^. The size of the sound horizon at this time is consequently imprinted as a 
bump in the two-point correlation function of the matter at late times. Assuming an FLRW evolution over the scale 
of the horizon at this time, this the proper size of the sound horizon at the drag epoch is approximately given by 
(assuming N^g = 3.04) 



F. Big Bang Nucleosynthesis and the Lithium problem 



5-fT [139]. 



G. The BAO 



ds = 
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Mpc 
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con- 
r\ from dif- 
Constraints 



from ^Li observations [139] in Galac- 
tic globular clusters and Galactic 
halo are shown separately, alongside 
"He [140] and ^He [141]. These agree 
with each other if ?7io ~ 4.5. On 
the other hand, D observations at 
high redshift (red) [142] and CMB re- 
quire r;io — 6. Bottom left we show 
how a varying radial profile for r;io 
(from ~ 4.5 at the centre to ~ 6 
asymptotically) can fit all the ob- 
servational constraints, for differing 
inhomogeneity scales. On the right 
are the nuclei abundances as a func- 
tion of z in an example model. This 
model may be considered an 'isocur- 
vature void' model. (From [75].) 



local measurements 




[Gpc] 



which is converted from [144] to make it purely local. In an FLRW model, this scale simply redshifts, and so can be 
used as a standard ruler at late times. In a void model, it shears into an axisymmetric ellipsoid through the differing 
expansion rates and H^. The proper lengths of the axes of this ellipse, when viewed from the centre, are given by 

^ll(^) = 4 = ,,^"^"1, y L^(z) = l^^'] .. = dAz)Seiz) , (29) 

a\\(td,r{z)) (l + z)iJ||(z) aj_[td,r[z)) 

where the redshift increment Sz{z) and angular size S9{z) are the corresponding observables. 
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FIG. 9: The BAO distance measure = {SO^Sz/ z^^'-^ shown 
compared to a ACDM model. The green and blue dashed lines 
are void models, whereas the purple and red lines are FLRW 
models. Clearly, the low redshift data favours FLRW over these 
models. (From [112].) 



Even without this shearing effect, the BAO can be used 
to constrain void models because it can be used to mea- 
sure H(z) through the 'volume distance': 



1/3 



(30) 



This quantity is given directly by current surveys. 

Thus, compared to from SNIa, the BAO pro- 

vide a complementary measurement of the geometry of 
the model, and in particular provide a probe of H{z). 
For the simplest types of voids with zero bang time 
and no isocurvature modes, the BAO are strongly in 
tension with the SNIa data [48, 49, 79, 82, 112] - see 
Fig. 9. Note that this assumes there are no compli- 
cations from the evolution of perturbations, and that 
scales evolve independently. While this is the case in 
FLRW, it is not in LTB where curvature gradients are 
present. Whether this is important to the analysis is 
yet to be shown (see below). 
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While the BAO arc indeed a restrictive test, it is clear that the constraints can easily be circumvented in the same 
way as the CMB. The bang time function can be used to free the constraint because it can be used to fix H\\{z) 
separately from dA{z) which is not the case if it is zero. Alternatively, we can use the freedom in = fb{i~) and 
V = v{^) to change dg as a function of radial shell about the observer [81, 112]. The BAO can then be interpreted 
as a measurement of these parameters in different shells around us. Similarly, radial changes in the primordial power 
spectrum can significantly affect these results [94]. While this might require some fine tuning to shift the BAO 
peak [112], it is not yet clear if this is a significant issue. 



H. Density perturbations 

An important open problem in inhomogeneous models is the modelling of striicture formation. This is important 
partly because it provides a means for distinguishing between FLRW and LTB. One example of where we might see 
an effect is in the peak in the two-point matter correlation function attributed to the Baryon Accoustic Oscillations 
(BAO). ft has been shown that if LTB perturbations evolve as in FLRW, then BAO can be decisive in ruling out 
certain types of voids [48, 49]. Whether this assumption is valid however requires a full analysis of perturbations. 

There have been three approaches so far: 

1. Using a covariant 1+1+2 formalism which was developed for gauge-invariant perturbations of spherically sym- 
metric spacetimes [145, 146]. The full master equations for LTB have not yet been derived, but some progress 
has been made in the 'silent' approximation, neglecting the magnetic part of the Weyl tensor [40, 72]. 

2. Using a 2+2 covariant formalism [147, 148], developed for stellar and black hole physics. The full master 
equations for LTB perturbations were presented in [53] (see also [149]). 

3. An N-body simulation has been used to study Newtonian perturbations of voids [90]. 

In FLRW cosmology, perturbations split into scalar, vector and tensor modes that decouple from each other, and 
so evolve independently (to first order). Such a split cannot usefully be performed in the same way in a spherically 
symmetric spacetime, as the background is no longer spatially homogeneous, and modes written in this way couple 
together. Instead, there exists a decoupling of the perturbations into two independent sectors, called 'polar' (or even) 
and 'axial' (or odd), which are analogous, but not equivalent, to scalar and vector modes in FLRW. These are based 
on how the perturbations transform on the sphere. Roughly speaking, polar modes are 'curl' free on S'^ while axial 
modes are divergence free. Further decomposition may be made into spherical harmonics, so all variables arc for a 
given spherical harmonic index £, and modes decouple for each £ - analogously to fc-modes evolving independently 
on an FLRW background. A full set of gauge-invariant variables were given by [148] who showed that there exists a 
natural gauge - the Regge- Wheeler gauge - in which all perturbation variables are gauge-invariant (rather like the 
longitudinal gauge in FLRW perturbation theory). Unfortunately, the interpretation of the gauge-invariant variables 
is not straightforward in a cosmological setting. 

Most of the interesting physics happens in the polar sector, so we will discuss that case, following [53]. The general 
form of polar perturbations of the metric can be written, in Regge- Wheeler gauge, as 

^d ^ 

= _ [1 + (2^ - X - >fi)Y] dt^ - -^^^dMr + [1 + (x + v)Y] , + alr^{l + ^F)d^2^ (31) 

VI - nr^ (1 - kH) 

where ri{t,r), x{t,r), <{>{t,r) and ?(t,r) are gauge-invariant variables. The notation here is such that a variable times 
the spherical harmonic Y has a sum over £,m, e.g., LpY ^ Y^^oYl^^^^Z^\^im{x^)Yim(.x°-), where a;" are coordinates 
on S"^, and = {t, r). The general form of polar matter perturbations in this gauge is given by 

p = p^rS(l + AF), (32) 

where v, w and A are gauge-invariant velocity and density perturbations and Hab is the metric perturbation in the 
part of the metric; a colon denotes covariant differentiation on the 2-sphere. The unit vectors in the time and 
radial directions are 

u^ = {l,0), = (^0,^^lE^y (33) 

The elegance of the Regge- Wheeler gauge is that the gauge-invariant metric perturbations are master variables for 
the problem, and obey a coupled system of PDEs which are decoupled from the matter perturbations. The matter 
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perturbation variables are then determined by the solution to this system. We outline what this system looks like for 
£ > 2; in this case ry = 0. The generalized equation for the gravitational potential is [53]: 



(34) 



The left hand side of this equation has exactly the form of the usual equation for a curved FLRW model, except that 
here the curvature, scale factor and Hubble rate depend on r. On the right, S^p is a source term which couples this 
potential to gravitational waves, x, and generalized vector modes, These latter modes in turn are sourced by (p: 



-X + X" -^H.x-2Wx' + 



X = S^{<;,ip), (35) 
4 + 2Huc; = -x'. (36) 



The prime is a radial derivative defined by X' = n^'V aX. 

The gravitational field is inherently dynamic even at the linear level, which is not the case for a dust FLRW 
model with only scalar perturbations. Structure may grow more slowly due to the dissipation of potential energy 
into gravitational radiation and rotational degrees of freedom. Since H± = H±{t,r), a± = a±{t,r) and k = K{r), 
perturbations in each shell about the centre grow at different rates, and it is because of this that the perturbations 
generate gravitational waves and vector modes. This leads to a very complicated set of coupled PDEs to solve for 
each harmonic £. 

In fact, things are even more complicated than they first seem. Since the scalar-vector-tensor decomposition does 
not exist in non-FLRW models, the interpretation of the gauge-invariant LTB perturbation variables is subtle. For 
example, when we take the FLRW limit we find that 



^ = -2* - 2'HV - 2^-^—^hr + ^h<^> + 



' + 1) - A{1 - Kr^) 
2^ 



(37) 



where 5* is the usual perturbation space potential, V is the radial part of the vector perturbation, and the ft-'s are 
invariant parts of the tensor part of the metric perturbation. Thus (p contains scalars, vectors and tensors. A similar 
expression for ? shows that it contains both vector and tensor degrees of freedom, while x ^ genuine gravitational 
wave mode, as may be seen from the characteristics of the equation it obeys. This mode mixing may be further seen 
in the gauge-invariant density perturbation which appears naturally in the formalism: 



87rGpA = -ip" -2Wip' + {H\\+2H^)ip + Wx' + HjlX 

(£-!)((, + 2) _ 



£{£ + 1) 



2 „2 



a I r 



+ 2HI+4HuHj_-8'kGp 



ix + v) 



2a\r'^ 



'X + 2HW + 2{Hu+H^)W^, 



where 



W = 



\/l — Kr^ 

a±r 



(38) 



(39) 



When evaluated in the FLRW limit the mode mixing becomes more obvious still: A contains both vector and tensor 
modes, while its scalar part is 



47rGaV A = V^* - SUdr^ - 3(7^^ _ 



(40) 



which gives the usual gauge invariant density fluctuation dp'"'' = 5p + drp{B — drE) [150]. Here, refers to the 
Laplacian acting on a 3-scalar. The fact that A is more complicated is because the gauge- invariant density perturbation 
includes metric degrees of freedom in its definition; gauge-invariant variables which are natural for spherical symmetry 
may not be natural for homogeneous backgrounds. A gauge-dependent A may be defined which reduces to Sp''^'' in 
the FLRW subcase, but its gauge-dependence will cause problems in the inhomogencous case. 

These equations have not yet been solved in full generality. We expect different structure growth in LTB models, 
but it is not clear what form the differences will take. It seems reasonable to expect that the coupling between scalars, 
vectors and tensors will lead to dissipation in the growth of large-scale structure where the curvature gradient is 
largest, as it is the curvature and density gradients that lead to mode coupling. In trying to use structure formation 
to compare FLRW to LTB models, some care must be taken over the primordial power spectrum and whatever early 
universe model is used to generate perturbations - since there is a degeneracy with the primordial power spectrum 
and the features in the matter power spectrum. 
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FIG. 10: Off-centre observers see 
anisotropy. 

Top: The main contribution to the 
total CMB anisotropy (a) is in the 
form of a dipole (b) with higher- 
order moments suppressed (c), (d). 
(From [151], for an LTB model with 
dark energy.) 

Bottom: There is also a dipole in the 
distance modulus, shown here for a 
large and small void at different red- 
shifts. (From [61].) 
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I. The Copernican problem: Constraints on the distance from the centre 



An ofF-centre observer will typically measure a large dipole in the CMB due to their peculiar velocity with respect 
to the CMB frame, which is non-perturbative in this context [4]. They will also measure a dipole in their local 
distance-redshift relation which is not due to a peculiar velocity effect nor a dipole in the Hubble law. Rather this 
first appears at 0{z^) in the distance-redshift relation through gradients in the expansion rate and the divergence of 
the shear - see Eq. (79) below. Combined constraints on the dipole in the SNIa data and the CMB would require us 
to be within ~ 0.5% of the scale radius of the void [164], without fine-tuning our peculiar velocity to cancel out some 
of the dipole. Others reach similar conclusions [4, 21, 25, 84]. 

An intriguing alternative view was presented in [79]. Although they reach the same conclusion as to how close to 
the centre the observer should be, they argue that if we're slightly off centre, then one would expect to see a significant 
bulk flow in the direction of the CMB dipole. Such a 'dark flow' has been tentatively measured [165, 166], and has 
not been accounted for within ACDM at present. 



J. Summary and interpretation of inhomogeneous models 

An inhomogeneous LTB void model, even if it over-simplifies nonlinear inhomogeneity at the background level, does 
produce some rather remarkable results. The apparent acceleration of the universe can be accounted for without dark 
energy, and the Lithium problem can be trivially solved. However, it seems that the simplest incarnation of such a 
model will not work: combining observables Hq, SNIa and the CMB reveals considerable tension with the data, with 
the main problem being a low Hq in the models [79]; add in (even rudimentary) kSZ and y-distortion constraints, and 
the situation is conclusive. Results from the BAO, though only indicative and not yet decisive (as they do not take 
into account structure formation on an inhomogeneous background), also signal considerable tension. 

Each of these observables point to non-trivial inhomogeneity at early times (or A of course) . Most models which 
are ruled out have the assumption of evolving from a homogeneous FLRW model. Primary observables provide much 
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weaker constraints if this restriction is removed, though it is still difficult to get a good fit using just the freedom of a 
bang time function [107]. But one can free up essentially any function which is assumed homogeneous in the standard 
model; in the context of inhomogeneous models, it doesn't make sense to keep them homogeneous unless we have 
a specific model in mind (perhaps derived from a model of inflation). Examples of such freedom include a radially 
varying bang time function, a radially varying primordial power spectrum (designed to have the required spectrum 
on 2-spheres perhaps), isocurvature degrees of freedom such as a varying baryon photon ratio or baryon fraction, and 
one can dream up more such as varying Neg. Indeed, taken at face value the lithium problem [139] can be interpreted 
as a direct measurement of an inhomogeneous isocurvature mode present at BBN in this context [75] this is actually 
the one observation which is potentially at odds with homogeneity. Most primordial numbers in the standard model 
are not understood well if at all, and if we remove slow roll inflation - as we must to make such a model in the first 
place - we remove significant motivation to keep them homogeneous. 

This suggests an important reverse-engineering way to handle such models. If we accept that presently any specific 
inhomogeneous model is essentially pulled from thin air, then we have to conclude that what we are really trying to 
do is to invert observables to constrain different properties of the model in different shells around us. Though it seems 
rather non-predictive, without an early universe model to create a void from it is really no different from making a 
map of the universe's history. This inverse-problem approach has been investigated in [45, 54, 80, 91, 122, 123, 125]. 
The idea is to specify observational data on our past lightcone, smooth it, and integrate into the interior. Whether 
this inverse problem is well- conditioned or not is crucial to the success of such an approach [110]. Nevertheless, this 
is a valuable strategy: If we specify data on our past lightcone and integrate into the interior, does it necessarily yield 
an FLRW model, or are there other solutions (perhaps without dark energy)? 

An important alternative view of these models is not to view them as an anti-Copernican alternative to dark 
energy, but rather to view them as the simplest probe of large scale inhomogeneity [106]. This is akin to considering 
Bianchi models as probes of large-scale anisotropy. We may therefore think of LTB models as models where we 
have smoothed all observables over the sky, thereby compressing all inhomogeneities into one or two radial degrees 
of freedom centred about us. In this interpretation, we avoid placing ourselves 'at the centre of the universe' in 
the standard way. Furthermore, constraints which arise by considering anisotropy around distant observers - the 
Goodman constraints - are perhaps removed from the equation; distant observers would see an isotropic universe too. 

In this sense, these models are a natural first step in developing a fully inhomogeneous description of the universe. 
There is a vital caveat to this interpretation, however: we must include dark energy in the model for it to be meaningful, 
which is almost never done (see [87, 151, 153] for a first attempt). If we do not, then we are implicitly assuming that 
consequences of the averaging, backrcaction and fitting problems really do lead to significant effects which solve the 
dark energy problem. That is, by averaging observables at some redshift over the sky we are averaging the geometry 
out to that redshift, which can have a non- trivial 'back-reaction' effect on our interpretation of the model [117]. This 
could conceivably look like dark energy in our distance calculations (perhaps even dynamically too [118]). If that 
were indeed the case, we could have a significant effective energy-momentum tensor which would be very different 
from dust, and it would not be simple to calculate observables as they would not necessarily be derivable from the 
metric. Hence, within this interpretation the dust LTB model would certainly be the wrong model to use. If one is 
to further peruse this idea, one might need to constrain deviations from homogeneity using the metric only, without 
resorting to the field equations at all (see [152] for further discussion). 

This is nevertheless in some ways the most natural way to place constraints on inhomogeneity. Yet, if large-scale 
inhomogeneity were present, we shall see in the next section that within GR it is challenging perhaps impossible - 
to reconcile it with the Copernican principle given the level of isotropy we observe. 

III. ROUTES TO HOMOGENEITY 

Considering a specific inhomogeneous solution to the EFE which violates the CP helps us consider the types of 

observables which can be used to demonstrate homogeneity. Ruling out classes of solutions as viable models helps test 
the Copernican principle provided we understand where they sit in the space of solutions, so to speak. Under what 
circumstances does the Copernican principle, combined with some observable, actually imply an FLRW geometry? 

In the case of perfect observables and idealised conditions quite a lot is known, as we discuss below. These results 
are non-perturbative, and do not start from FLRW to show consistency with it. For the case of realistic observables 
in a lumpy universe, however, details are rather sketchy, with only one case properly considered. 

Many of these results rely on the following theorem [154]: 

The FLRW models: For a perfect fluid solution to the Einstein field equations where the velocity field of the fluid 
is geodesic, then the spacetime is FLRW if either: 

- the velocity field of the source is shear-free and irrotational; or. 
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- the spacetime is conformally flat (i.e., the Weyl tensor vanishes). 

The 'perfect fluid source' here refers to the total matter content and not to the individual components, and is 
necessarily barotropic if either of the conditions are met. So, for example, the matter could be comoving dark matter 
and baryons, and dark energy in the form of a scalar fleld with its gradient parallel to the velocity of the matter. 

A. Isotropy of the CMB 

What can we say if the CMB is exactly isotropic for 'fundamental observers'? This is the canonical expected 
observable which intuitively should imply FLRW with the CP. It does, usually, but requires assumptions about the 
theory of gravity and types of matter present. The pioneering result is due to Ehlers, Geren and Sachs (1968) [155]. 
Without other assumptions we have: 

_|_ 

[EGS] Radiation isotropy > conformally stationary: In a region, if observers on an expanding congruence 

u° measure a collisionless radiation field which is isotropic, then the congruence is shear-free, and the expansion 
is related to the acceleration via a potential Q = — jlnp^: = D^Q, O = SQ; the spacetime must be 
conformal to a stationary spacetime in that region. 

That doesn't tell us a great deal, but including geodesic observers changes things considerably. The original EGS 
work assumed that the only source of the gravitational field was the radiation, i.e., they neglected matter (and they 
had A = 0). This has been generalised over the years to include self-gravitating matter and dark energy [78, 156-161], 
as well as for scalar-tensor theories of gravity [162] : 

[EGS-H] Radiation isotropy with dust > FLRW: In a region, if dust observers on an expanding congru- 
ence u°' measure a collisionless radiation field which has vanishing dipole, quadrupole and octopole, and non- 
interacting dark energy is the form of A, quintessence, a perfect fiuid or is the result of a scalar-tensor extension 
of GR, then the spacetime is FLRW in that region. 

The dust observers are necessarily geodesic and expanding: 

Aa^O, e>0. (41) 

Because the dust observers see the radiation energy flux to vanish (the dipole), is the frame of the radiation 
also. The photon distribution function f{x,p, E) in momentum space depends on components of the 4-momentum p" 
along u"-, i.e., on the photon energy E = —UaP"', and, in general, the direction e", and may be written in a spherical 
harmonic expansion as (see the appendix) 

oo 

/ = ^F^,e^^ (42) 
e=o 

where the spherical harmonic coefficients Fai are symmetric, trace- free tensors orthogonal to u", and Ai stands for 
the index string aia2 ■ ■ ■ ag. (In this notation, e^^''' are a representation of the spherical harmonic functions.) The 
dust observers measure the first three moments of this to be zero which means 

Fa = Fab = Fabc = 0. (43) 

In particular, as follows from Eq. (131), the momentum density (from the dipole), anisotropic stress (from the 
quadrupole), and the radiation brightness octopole vanish: 

g« = 7rf = n"*"^ = . (44) 

These are source terms in the anisotropic stress evolution equation, which is the ^ = 2 case of Eq. (137). In general 
fully nonlinear form, the -k^" evolution equation is 

o lo 

+ ^<7e<V^>^ + ^D,n«^^ - ^^a^W"^ = 0. (45) 
7 oo oio 



18 



Eq. (44) removes all terms on the left except the third and the last: 

(21p,/i„V-47rn„6^'^)aed = 0, (46) 

which imphes, since Habcd is trace-free and the first term consists of traces, shear-free expansion of the fundamental 
congruence: 

<Tab = 0. (47) 

We can also show that is irrotational as follows. Together with Eq. (41), momentum conservation for radiation, 
i.e., Eq. (126) with I = r, reduces to 

DaPr = 0. (48) 

Thus the radiation density is homogeneous relative to fundamental observers. Now we invoke the exact nonlinear 
identity for the covariant curl of the gradient, Eq. (104): 

CmlDaPr = -"^PrOJa 0/9^^0 = 0, (49) 

where we have used the energy conservation equation (125) for radiation. By assumption 6 > 0, and hence we deduce 
that the vorticity must vanish: 

OJa = 0- (50) 

Then we see from the curl shear constraint equation (118) that the magnetic Wcyl tensor must vanish: 

Hab = 0. (51) 

Furthermore, Eq. (48) actually tells us that the expansion must also be homogeneous. From the radiation energy 
conservation equation (125), and using Eq. (44), we have 9 = —3pr/'ipr. On taking a covariant spatial gradient and 
using the commutation relation Eq. (105), we find 

D„e = 0. (52) 

Then the shear divergence constraint, Eq. (117), enforces the vanishing of the total momentum density in the matter 
frame, 

/ I 

The second equality follows from Eq. (AID) in [163], using the fact that the baryons. CDM and dark energy (in the 
form of quintessence or a perfect fluid) have vanishing momentum density and anisotropic stress in their own frames, 
i.e., 

qr = = 7rr\ (54) 

where the asterisk denotes the intrinsic quantity (see Appendix A). If we include other species, such as neutrinos, 
then the same assumption applies to them. Except in artificial situations, it follows from Eq. (53) that 

v'} = 0, (55) 

i.e., the bulk peculiar velocities of matter and dark energy [and any other self-gravitating species satisfying Eq. (54)] 
are forced to vanish - all species must be comoving with the radiation. 

The comoving condition (55) then imposes the vanishing of the total anisotropic stress in the matter frame: 

^ab ^J2^f ^ J2^j{p*+p*)v\"vf =0, (56) 

where we used Eqs. (All) in [163], (54) and (55). Then the shear evolution equation (113) leads to a vanishing electric 
Weyl tensor 

Eab = 0. (57) 
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Equations (53) and (56), now lead via the total momentum conservation equation (111) and the S-divergence con- 
straint (119), to homogeneous total density and pressure: 

D„p = = D„p. (58) 

Equations (41), (47), (51), (52), (53), (56) and (58) constitute a covariant characterisation of an FLRW spacetime. 
This establishes the EGS result, generalised from the original to include self-gravitating matter and dark energy. It 
is straightforward to include other species such as neutrinos. 

The critical assumption needed for all species is the vanishing of the intrinsic momentum density and anisotropic 
stress, i.e., Eq. (54). Equivalently, the energy- momentum tensor for the /-component should have perfect fluid form 
in the /-frame (we rule out a special case that allows total anisotropic stress [161]). The isotropy of the radiation 
and the geodesic nature of its 4-velocity - which follows from the assumption of geodesic observers - then enforce the 
vanishing of (bulk) peculiar velocities Vj. Note that one does not need to assume that the other species are comoving 
with the radiation - it follows from the assumptions on the radiation. A similar proof can be used for scalar-tensor 
theories of gravity, although it is somewhat more involved [162]. 

It is worth noting that we do have to assume the the dark matter and dark energy are non-interacting. If we do 
not, we cannot enforce the radiation congruence to be geodesic because the observers may not be, and one is actually 
left only with a very weak condition: only that the spacetime is conformally stationary [155, 160, 167]. 

In summary, the EGS theorems, suitably generalised to include baryons and CDM and dark energy, are the most 
powerful basis that we have - within the framework of the Copernican Principle - for background spatial homogeneity 
and thus an FLRW background model. Although this result applies only to the 'background Universe', its proof 
nevertheless requires a fully nonperturbative analysis. 

B. Blackbody spectrum of the CMB 

The EGS results rely only on the isotropy of the radiation field, and do not utilise its spectrum. Of course, no 
geometry can affect the spectrum of the CMB because the gravitational influence on photons is frequency independent 
(except at very high frequencies). However, the fact that the CMB is a nearly perfect blackbody tells us much about 
the spacetime when there are scattering events present. The Sunyaev-Zel'dovich (SZ) effect is due to the scattering 
of CMB photons by charged matter, and has already been shown to be a powerful tool for constraining radial 
inhomogeneity within the class of cosmological models constructed from the LTB solutions, as discussed in Sec. HE. 
It has recently been shown [168] that under idealised circumstances similar to the EGS theorems above the SZ effect 
can actually be used as a proof of FLRW geometry for one observer without requiring the Copernican principle at 
all, thus extending Goodman's tests to arbitrary spacetimes [135]. 

[GCCB] Isotropic blackbody CMB > FLRW: An observer who sees an isotropic blackbody CMB in a uni- 
verse with scattering events can deduce the universe is FLRW if either double scattering is present or they 
can observe the CMB for an extended period of time, assuming the CMB is emitted as a blackbody, and the 
conditions of the EGS theorem hold. 

The observational effect of the scattering of CMB photons by baryonic matter is usually referred to in the literature 
as the SZ effect [136, 137], and is often divided into two different contributions; the thermal SZ effect (tSZ) [136] 
and the kinematic SZ effect (kSZ) [137]. The kSZ effect causes a distortion in the spectrum of the reflected light 
due to the anisotropy seen in the CMB sky of the scatterer, and maintains the same distribution function it had 
before the scattering event (all other changes being encapsulated in the tSZ). For the case of blackbody radiation this 
corresponds solely to a change in temperature of the scattered radiation. 

Thus, an observer who sees an exactly isotropic CMB in a universe where scattering takes place, can deduce that 
the scatterers themselves must also see an isotropic CMB, provided that decoupling emits the CMB radiation as an 
exact blackbody. The proof of this relics on the fact that blackbody spectra of differing temperatures cannot be added 
together to give another blackbody. This mechanism therefore provides, in effect, a set of mirrors that allows us to 
view the CMB from different locations [135]. 

Is this enough to deduce FLRW geometry? Not on its own. Such an observation gives a single null surface on 
which observers see an isotropic CMB. This allows us to use a much weaker version of the Copernican principle than 
used in the EGS theorems (which assume isotropy for all observers in a 4-dimensional patch of spacetime, and deduce 
homogeneity only in that patch) to deduce homogeneity. For example, if all observers in a spacelike region see an 
isotropic blackbody CMB when scattering is present, then the spacetime must be homogeneous in the causal past of 
that region. 
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FIG. 11: The SZ effect provides 
information about the isotropy 
of the CMB sky at other points 
on our past hghtcone. It can 
also provide us with information 
about parts of the last scatter- 
ing surface that would otherwise 
be inaccessible to us. Multiple 
scattering events provide further 
information about the CMB sky 
at other points within our causal 
past. (From [168].) 
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The lone observer can say more, however [168]. If there are scattering events taking place throughout the universe, 
then each 'primary scatterer' our observer sees must also be able to deduce that scatterers on their past nuUcone see 
an isotropic CMB, or else their observations of a blackbody spectrum would be distorted. Consequently, the spacetime 
must be filled with observers seeing an isotropic CMB, and one can then use the ECS theorems to deduce FLRW 
geometry. Hence, under highly idealised conditions, a single observer at a single instant can deduce FLRW geometry 
within their past lightcone. Alternatively, the observer could wait for an extended period until their past nuUcone 
sweeps out a 4-D region of spacetime. If no kSZ effect is measured at any time, then they can infer that their entire 
causal past must also be FLRW. 



C. Local observations 



Instead of using the CMB, local isotropy of observations can also provide a route to homogeneity. Adopting the 
Copernican Principle, it follows from the 'lightcone-isotropy implies spatial isotropy' theorem that if all observers see 
isotropy then the spacetime is isotropic about all galactic worldlines - and hence spacetime is FLRW. 

+ CP 

Matter lightcone-isotopy > FLRW: In an expanding dust region with A, if all fundamental observers mea- 
sure isotropic area distances, number counts, bulk peculiar velocities, and lensing, then the spacetime is FLRW. 

In essence, this is the Cosmological Principle, but derived from observed isotropy and not from assumed spatial 
isotropy. Note the significant number of observable quantities required. Using the CP, we can actually give a much 
stronger statement than this, based only on distance data. An important though under-recognised theorem due to 
Hasse and Perlick tells us that [169] 

+ CP 

[HP] Isotropic distances > FLRW: If all fundamental observers in an expanding spacetime region measure 

isotropic area distances up to third-order in a redshift series expansion, then the spacetime is FLRW in that 
region. 

The proof of this relies on performing a series expansion of the distance-redshift relation in a general spacetime, 
using the method of Kristian and Sachs [170], and looking at the spherical harmonic multipoles order by order. We 
illustrate the proof in the 1-1-3 covariant approach (in the case of zero vorticity). Performing a series expansion in 
redshift of the distance modulus, we have, in the notation of [171], 

m-M -25 = 51ogioz-51ogio X"if''VaU6|„ + ^logioe 

_ " RabK'^K'' _ 3(j^°i^''j^°VaVbMe)' K'' K'> K^'V cUd 

where 

^""^^rr K-\o^ -u^ + e''\o^ (60) 



7? ^ -|-0(z3),(59) 
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denotes a past-pointing null vector at the observer O in the direction e°. When fully decomposed into their projected, 
symmetric and trace-free parts, products of e's represent a spherical harmonic expansion. Thus, this expression views 
the distance modulus as a function of redshift on the sky, with a particular spherical harmonic expansion on a sphere 
of constant redshift. (The inverse of this expression has coefficients which have a spherical harmonic interpretation 
on a sphere of constant magnitude.) Comparing with the standard FLRW series expansion evaluated today, we define 
an observational Hubble rate and deceleration parameter as 



H 



obs I 



^obs I 



^ 10 {K^RPVaUp)^ 
We can also give an eff'ective observed cosmological constant parameter from the 0(-2^) term: 



(1 - gg*^^) - 5 



12(_H-obs)2 



6(iJ°bs)3 



(61) 

(62) 
(63) 



The argument of [169] relies on proving that if all observers measure these three quantities to be isotropic then the 
spacetime is necessarily FLRW. In a general spacetime [171] 



tjobs 
-"0 



a J) 



(64) 



where A^e" is a dipole and aabe-°'e^ is a quadrupole. Hence, if all observers measure H^^^ to be isotropic, then 
Cab = = ^a- In a spacetime with isotropic Hq^^ the generalised deceleration parameter, defined on a sphere of 
constant redshift, is given by [172, 173] 



2qa - 3V„e 



(65) 



If the dipole of this term vanishes then we see from Eq. (117) that the energy flux must vanish as well as spatial 
gradients of the expansion. Excluding models with unphysical anisotropic pressure, Eq. (114) then shows that the 
electric Weyl tensor must vanish, and it follows that the spacetime must be FLRW. The more general proof in [169] 
uses flf^^ to show that the vorticity must necessarily vanish along with the anisotropic pressure. 



D. The Hubble rate on the past lightcone 



Measurements of the Hubble rate on our past lightcone can provide an important route to homogeneity assuming 
the Copernican principle. In a general spacetime, the spatial expansion tensor may be written as 

©ab = Qhab + (Tab ■ (66) 

However, we do not measure exactly this as our observations are made on null cones. An observer can measure the 
expansion at a given redshift in three separate directions: radially, and in two orthogonal directions in the screen 
space. The observed radial Hubble rate at any point is a generalisation of our observed Hubble constant above, 

ffjl (z; e) = K'^K^'VaUb = ^6 - A^e" + (7abe"e'' . (67) 
The Hubble rates orthogonal to this may be found by projecting into the screen space, Nab = 9ab—KaKb—KaUb—UaKb: 

H±ab{z; e) = ^Na^Mb'^VeUd = - la,,e'e^^ Nab + (e^o;'^) + \ [Na'Nb'' - h^ab^"^^ cjab ■ (68) 

The trace of this gives the areal Hubble rate 

H^{z; e) = H^abN'''' = ^6 - ^a„6e«e'' , (69) 

which implies the spatial volume expansion rate is given via the observable expansion rates as: 

e{z;e)=H\\{z;e) + 2H^{z;e) + Aae''. (70) 
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If we measure H^^i, to be rotationally symmetric on the celestial sphere in each direction e° and at each z, and 
(z; e) = H±{z; e) then this, on its own, is not enough to set the shear, acceleration and rotation to zero, even on 
our past lightcone. However, we can see how measuring the Hubble rate can lead to homogeneity, when we apply the 
Copernican principle: 

Isotropic Hubble rate > FLRW If all observers in a perfect fluid spacetime measure: 

- H±ab to be rotationally symmetric on the screen space at each z ; and, 

- H\\{z;e) = Hj_{z]e) , 

then the spacetime is FLRW. 

The proof of this is straightforward: the first condition implies that = aab = and the second that A"- = 0, which 
implies FLRW from the theorem above. 

This forms the basis of the Alcock-Paczynski [174] test in a general spacetime. An object which is spherical at one 
instant will physically deform according to Eq. (66), if it is comoving (such as the BAO scale). A further efl'ect is that 
it will be observed to be ellipsoidal with one of the principle axes scaled by H^^ and the other two by the inverse of 
the eigenvalues of H±ab- 



E. Ages: absolute and relative 

Measurement of the ages of objects - both absolute and relative - on our past lightcone provides important 
information which can be used to deduce homogeneity. 

Neighbouring lines on the matter congruence Ua = dot, measuring proper time t, may be thought of as connected 
by null curves fc° (which are past pointing in our notation above). An increment of redshift on a null curve is related 
to an increment of proper time on the matter worldlines as 

^ = I (71) 

dz (l + z)i/||(z;e) ^ ' 

This gives the age difference of objects observed on the past lightcone, in a redshift increment dz. Over a finite 
redshift range, an observer at the origin may determine the age difference between objects A and B as 



tA — ts 



r ^ (72) 

A, (l + z)if||(z;e)' ^ > 



where the integral is along the null curve connecting A and B. In a general spacetime, the absolute age of an object 
is given by the time interval from the big bang, which may not be homogeneous. Therefore, 



ft 

T 



/ dt = t-tBB, (73) 



where the integral is along the worldline of the object. So, in terms of absolute ages we have 

+ z)Hv.{z;e) 



rZB 

TA-TB= / , , ^ - tBB{A) + tBB{B) . (74) 

JzA U- - 



Measurements of ages, then, provide two important insights into inhomogeneity: firstly they probe the radial 
Hubble rate; secondly they give a direct measure of the bang time function, which arrises because surfaces of constant 
time may not be orthogonal to the matter. In a realistic model, this could represent the time at which a certain 
temperature was attained (so switching on a given type of cosmological clock). More generally, the 'big bang' need 
not be homogenous, and ages provide a mechanism to probe this, given a separate measurement of Hu [175]. 



F. Does 'almost' isotropy imply 'almost' homogeneity? 

The results above arc highly idealised because they assume perfect observations, and observablcs such as isotropy 
which are not actually representative of the real universe. In reality, the CMB temperature anisotropies, though small 
are not zero; local observations are nearly isotropic, but not to the same degree as the CMB. Are the above results 
stable to these perturbations? 

The key argument is known as the almost-EGS theorem [158, 176-179]: 
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CP 



> almost FLRW: In a region of an expanding Universe with dust and cos- 



[SME] Almost isotropic CMB 

mological constant, if all dust observers measure an almost isotropic distribution of collisionless radiation, then 
the region is almost FLRW, provided certain constraints on the derivatives of the multipoles are satisfied. 

The starting point for this proof is to assume that the multipoles of the radiation are much smaller than the 
monopole, which is exactly what we measure: 



pr 



for 



1,...,4 



(75) 



where ©(e^) is a smallness measure. However, we can see from Eq. (45) that we need further assumptions on the 
derivatives of the multipoles to prove almost homogeneity: 



in„ 



= 0{er) for £ = !,..., 3 



(76) 



The proof proceeds as in the exact case above, but with = replaced by = 0{er) (except for Ag, = exactly because 
the observers are dust). An 'almost FLRW condition is then arrived at: 



kinematics : \aab\ = |wa| = |Dc(Tab| = IDbWa] = |Da6| = • 
curvature : \Eab\ = \Hab\ = IDc^^afcl = \'DcHab\ = 0{er) , 



0{er) 



(77) 
(78) 



in the region where the CMB is close to isotropic. Outside that region the spacetime need not be close to FLRW [180]. 

This proof is an important attempt to realistically arrive at a perturbed FLRW model using the CP and observables. 
It has been criticised as a 'realistic' basis for ncar-homogcncity due to the reliance of the assumptions on the spatial 
gradients of the multipoles, Eq. (76) [181]-^. However, it seems reasonable that the gradients of the very low multipoles 
are small compared to their amplitude (i.e., it seems reasonable that the CMB power spectrum does not change a lot 
as we move from observer to observer). This is because the multipoles peak in power around the scale they probe, 
which for the low multipoles required for the theorem are of order the Hubble scale. So, while they may change rapidly 
on small scales, the power of such modes is significantly diminished. Whether or not such criticisms are justified, the 
fact we have to make such assumptions means that the observational case in the real universe is much harder than 
the exact results imply. How could we observe spatial gradients of the octopole of the CMB? 

One route around this is to combine local distance information as well, and try for an almost- HP theorem. Consider 
the case of irrotational dust. Let us assume that all observers measure an approximately isotropic distance-redshift 
relation, and that the multipoles are bound by 0{ed). Clearly, we have, if this is the bound for all observers, 
Ca6 = 0{ed). Using the fully non-linear expression for the 0(-2^) term requires near-isotropy of 



K-K^K-VaVbUo =\{p + ip) + -\a + CTabo''^ + 
6 3 3 



Eab + gQcTab + ^G^aUbc 



1 2 

-Dae+ -dlVCTa 
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(79) 



which, together with Eq. (116), and assuming that time derivatives of 0(ed) quantities are 0(ed), give the conditions 



I dive 



ID 



IDaOl = \Egb\ = IdivKl = |divff„| = 0(ed) . 



(80) 



Yet, we still cannot arrive at the almost-FLRW condition as there is not enough information to switch off 'curl' degrees 
of freedom: i?a6) curl Sab, and curl (Tab are unconstrained. Polarisation of the CMB could be used to constrain these, 
but it would be interesting to see how everything fits together. 

A critical issue with these arguments lies in the question of what we mean by 'almost-FLRW, and what '= 0(e)' 
means. A sensible definition might be a dust solution which is 'almost' conformally flat. (Exact conformal flatness 
implies exactly FLRW.) Or, as used in the almost-EGS case, a set of conditions on all 1+3 irreducible vectors and 
PSTF tensors (which are necessarily gauge-invariant in the covariant approach) having small magnitudes: 



IX, 



ah-- 



^ Xg^...,X-^-- = 0(e). 



(81) 



^ Note that [181] is discussing something slightly different to the almost EGS theorem. In almost EGS, the assumption is that all observers 
measure nearly isotropic radiation; in [181] the assumption is that p'^\'aE is almost isotropic, but this is not the same thing. In the exact 
case this condition implies that the acceleration must be zero independent of the matter — this is not necessary for isotropic radiation, 
cf the radiation isotropy condition above. 
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Most quantities at first order average to zero in a standard formulation, so taking magnitudes is important. One 
of the conditions which does not rely on coordinates might therefore be EabE""^ = 0{e^) in the notation above. 
Unfortunately, this is non-trivial. Consider a linearly perturbed flat ACDM model in the Poisson gauge. Then, 



Evaluating the ensemble average of this gives its expectation value, which, assuming scale-invariant initial conditions 
gives [173] 



Here, fcuv is the wavenumber of some UV cutoff, and fceq is the wavenumber of the equality scale. The ensemble 
average has been evaluated assuming ergodicity via a spatial integral over a super-Hubble domain. We have a scalar 
which is 0(1) times a term which diverges in the UV. The divergence here represents modes which arc smoothed over, 
and seems to be a necessary condition for writing down an approximately FLRW model. Eq. (83) is certainly not 
0{e) in the normal sense of the meaning. This is because, roughly speaking, the Weyl tensor has a large variance even 
though we would like it to be 'small' for our covariant characterisation of FLRW. This is true also for most covariant 
objects with indices - their magnitudes are actually quite large in a perturbed FLRW model! Products of quantities 
such as <Jab<y°'^ and E^i^g^a^)'^ act as sources in the 1+3 equations. 

Instead we might try to define almost FLRW in the conventional perturbative sense. That is, we write the metric 
in the Newtonian gauge and require that the potential its first derivatives, and the peculiar velocity between the 
matter frame and the coordinate frame is small. This is often claimed to be a sufficient condition for a spacetime to 
be close to FLRW [182]. While this is fine in certain contexts, it is not necessarily so for large-scale inhomogeneities 
we are concerned with here. The LTB models provide a counter-example, as they can also be written as perturbed 
FLRW [183], yet arc clearly not 'almost FLRW' in the sense wc arc interested in. 

Consequently, a robust, non-perturbative, covariant definition of almost-FLRW is lacking and we are forced into 
the realm of the averaging problem: to define almost-FLRW, we must smooth away power on scales smaller than a 
few Mpc. How should this be done covariantly and what are the implications [117, 118]? 

IV. NULL HYPOTHESES FOR FLRW AND TESTS FOR THE COPERNICAN PRINCIPLE 

As we have seen, it is rather difficult to conclude homogeneity even given ideal observations. Indeed, it is subtle 
to robustly deduce (approximate) spherical symmetry about us given (near) isotropy of observations. It is rather 
surprising how sparse our current observations are in comparison to what is required from the theorems above, one 
of which includes observing transverse proper motions for an extend period! 

Nevertheless, the results above allow us to formulate some generic tests of the Copernican and cosmological prin- 
ciples, and the underling assumption of an FLRW geometry and its accompanying observational relationships - i.e., 
the usual observational relationships which arc derived assuming there are no issues from averaging. Wc refer to this 
as the 'FLRW assumption' below. (If there are non-trivial effects associated with averaging, then some of the tests 
below can also be used to signal it - see e.g., [184] for some specific examples.) 

The types of tests we consider range in power from focussed tests for the concordance model (flat ACDM) , to generic 
null tests which can signal if something is wrong with the FLRW assumption under a wide variety of circumstances 
irrespective of dark energy, initial conditions or theory of gravity. In this sense, they formulate our understanding of 
the Copernican principle as a null hypothesis which can be refuted but not proven. One of the conceptually important 
issues with some of these tests is that they can, in principle, be utilised by simply smoothing observed data, and so 
do not require an underlying model to be specified at all. These tests should be considered as additional to checking 
various observables for isotropy, which we assume below. 



(82) 




(83) 



A. 



Tests for the concordance model 



A simple consistency test for flat ACDM may be formulated easily. In a flat ACDM model we can rearrange the 
Friedmann equation to read [185, 186] 
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FIG. 12: Left: ffm{z) obtained us- 
ing distances for different FLRW 
models. The figure shows the 
range of behaviour from curvature 
in each fan of curves (key, top left). 
The three grey fans show how dif- 
fering r2„i values interact with cur- 
vature for ACDM (key, top right). 
The effect of changing w by a con- 
stant is illustrated in the red and 
brown fans. (From [187].) 
Right: Constraints on ffm{z) 
(called Om{z) in [186]) from 
the latest Wigglez BAO data. 
(From [188].) 
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where h{z) = H{z)/Hq is the dimensionless Hubble rate. Viewed in terms of the observable functions on the rhs, 
these equations tells us how to measure ftm from H{z) or D'{z). Furthermore, if flat ACDM is correct the answer 
should come out to be the same irrespective of the redshift of measurement. This provides a simple consistency test 
of the standard paradigm - deviations from a constant of ^m(^) at any z indicates either deviations from flatness, or 
from A, or from the FLRW models themselves [185, 186]. See Fig. 12. 

To implement it as a 'consistency test' one firsts smooths some data in an appropriate mo del- independent way, and 
then constructs the function ffm{z) to determine if it is consistent with the null hypothesis that it should be constant. 
See [187] for a discussion of how to do this with SNIa data, and [188] for an example using the BAO to measure H{z). 

More generally, if we don't restrict ourselves to 0,^ — 0, we have that [189]'^ 

nk = T(z){2{l-[l + zf)D" + -W{D'^ -l){l + zf] = d^\z) ^^°^> const. (85) 
= 2T(z) { [(1 + zf -D""-!] D" - [D'^ - l) [(1 + z)D' - D]] = fffX^) const. (86) 

where 

T(z)-i ^ 2 [1 - (1 + zf\ D'^D" - {(1 + z) [(1 + zf - 3(1 + z) + 2] i?'^ - 2 [l - (1 + zf] DD' - 3(1 + zfo'^] D' . 

The numerator of the formula for Vl^ forms the basis of the test presented in [185]. Again, these formulae for and 
fifc can be used to test consistency of the ACDM model. Note that each of these tests depend only on D{z), and not 
on any parameters of the model. 



B. Tests for FLRW geometry 

1. Hubble rate(s) and ages 

In an FLRW model, the two Hubble parameters we have discussed, i/y (z; e) and H±{z; e) measured in any direction 
must be the same at the same redshift. That is, 

{z;e) = Hj_ (z; e) = H{z) . (87) 

In a typical LTB model, for example, this simple relation is violated (with the exception of a fine-tuned sub-class of 
models of course), and so this forms a basic check of the FLRW assumption, and Copernican principle. Thus we have 
an isotropic expansion test: 

Jf{z) = H\\{z;e)-H^{z-e) 0. (88) 



^ Thanks to Marina Seikel and Sahba Yahya for correcting an error in these formula in [189]. 
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An important question here is how H±^ can be measured. In [38, 62, 72, 88, 189] it was shown that, in LTB models 
for a central observer, the 'redshift drift' - the change in redshift of a source measured over an extended time - is 
determined by the angular Hubble rate: 

z{z) = {l + z)Ho-H^{z). (89) 

Although it is a rather sensitive observable [190], it gives an important consistency check for the standard model. 
Another possibility to measure Hj_{z) is the polarised SZ effect (see below). 

Measuring ages of objects is notoriously difficult, yet can provide important insights as it probes the geometry 
inside our past null cone. At the simplest level, the relative ages of objects can be measured between redshift bins 
giving the Hubble rate. At a redshift z, the Hubble rate is may be found from Eq. (71) if we know the age difference 
between two nearby objects, dt separated by 5z: 

The fossil record of galaxies provides, in principle, the absolute ages of galaxies as a function of z assuming a model 
of galax;y evolution. In an FLRW model, this constrains the line of sight integral of the Hubble rate as [175] 

/•ZgBlaxy 2 FLRW 

y{z) = r.. - r,.u.y - I il + ,)H,iz) ' ° " ^''^ 

This has yet to be implemented. 

The dipole in the distance data can measure H{z) in a perturbed FLRW context [191]. It would be interesting to 
see in a more general context whether this is Hu (z), H±{z) or Q{z), thereby determining how it can be used for Jf{z). 



2. Curvature test 



A rather general test of FLRW geometry may be found by considering the distance-redshift relation in an open or 
closed FLRW model. Consider the dimensionless comoving distance: 

D{z) = (1 + zWAiz) = sin [y^kf^ ^) , (92) 

where h{z) = H{z)/Ho and cIa is the area or angular diameter distance. Rearranging Eq. (92), we have the curvature 
parameter today given by, for any curvature [37, 192] 

[h{z)D'{z)f - 1 ,^ . s FLRV^f ^ 

= \d{z)]2 = ^kiz) > const. (93) 

On the face of it, this gives a way to measure the curvature parameter today by combining distance data with Hubble 
rate data, irrespective of the redshift of measurement. In FLRW this will be constant as a function of z, independently 
of the dark energy model, or theory of gravity. Alternatively, we may re-write this as the condition that [37] 

^{z) = l + h''{DD" -D'^)+hh'DD' 0, (94) 

by virtue of Eq. (93). In more general spacetimes this will not be the case. In particular, in LTB models, even for a 

for a central observer, we have "^(z) 7^ 0, or ^k{z) 7^ const. 

This tells us that in all FLRW models there exists a precise relationship between the Hubble rate and distance 
measurements as we look down our past null cone. This relationship can be tested experimentally without specifying 
a model at all, if we reconstruct the functions H{z) and D{z) in a model independent way and independently of each 
other. This provides a model-independent method by which to experimentally verify the Copernican assumption, and 
so verify the basis of the FLRW models themselves (see Fig. 13). 

The curvature test Eq. (93) has been implemented in [187] using a mixture of SNIa, ages [to give H{z)], and BAO 
data; constraints are weak at present (Fig. 13). This is because the parameter is unstable at low redshift, while at 
high-z the SNIa data is very sparse. A variant of this was presented in [193] where it is reformulated as an integral test 
to measure curvature with much better results. Consistency between D{z) and H{z) measurements was explored and 
found in [194]. The dipole in the distance data can measure H{z) [191] so can be used directly in the curvature tests 
'^{z) and ^k{z)- This means that these curvature tests can be formulated from distance data alone if we interpret 
D{z) as the monopole of the distance-redshift relation in Eq (94), though this has not yet been investigated. 



27 



Reconstructed Q,, 



Reconstructed Q 




FIG. 13: Left: Consistency functions, shown for anti-Copernican models. The functions ^rn{z), ^k{z), -^{z) = 
[numerator of ^i^\z)] [185] and "^{z) are shown for a variety of void models. Deviations from constant indicate inconsis- 
tency with the standard paradigm, but only i?k{z) 7^ const, and '^(2) 7^ definitely indicates the FLRW assumption is wrong. 
(From [60].) 

Right: Constraints on &k{z) from a combination of SNIa and H(z) data from the relative ages of LRGs. This function is most 
usefully constrained at high redshift (blow-up). (From [187].) 



CMB 



Goodman's tests allow us to directly 'observe' the CMB from other locations. Backed up with the EGS theorems, 
and the scattering theorem discussed in Sec. IIIB, they help us to deduce homogeneity. Different multipoles can be 
probed in different ways, and so provide independent consistency conditions on the geometry at late times. If we 
consider fundamental observers comoving with the matter flow (i.e., without boosting to the CMB frame) we have: 

Monopole In FLRW, the monopole temperature around distant observers redshifts as (1 -I- z), which is not the case 
in an inhomogeneous geometry. (The direct temperature along a line of sight redshifts like this of course, but 
the monopole does not.) Thus we have, assuming standard physics. 



{Tcmb){z) 
To{l + z) 



- 1 



FLRW 



> 0, 



(95) 



where (T'cmb)('2) is the monopole temperature observed at redshift z in some frame w", and Tq is our CMB 
temperature. This can be measured, for example, via the thermal SZ effect [135] and via atomic fine-structure 
transitions in cool absorption-line systems along the line of sight to high-redshift quasars [195]. Constraints 
are reasonably tight already: [196] find that possible deviations of the form {Tcmb){z) = To(l -I- z)^~^ to be 
/3 = 0.004 ± 0.016 up to z ^ 3 from SZ measurements. 

Dipole Measuring the radial dipole around other observers using the kSZ effect constrains a mixture of their local 
peculiar velocity (which has a component to the mean matter frame as well as to the CMB frame) as well as 
spatial temperature gradients in the monopole (and higher-order effects) along the line of sight (cf Eq. 126 with 
I = r). Obviously this must be small in the standard model, so a strong kSZ signal would violate homogeneity. 
The integrated kSZ signal is consistent with the standard model [86] (see Sec. HE), but there have been 
measurements of very large dipoles measured around X-ray clusters [166], which has been interpreted as a 'dark 
flow', and is not consistent with ACDM. 

The y-distortion also depends on the total dipole which, if other contributions to the y-distortion can be 
accounted for, offers the possibility of inferring transverse velocities by comparing it with the kSZ signal (which 
depends only on the radial dipole). Transverse velocities are a key ingredient for determining isotropy of the 
lightcone, as well as for determining H±{z), so constraining these are a key test for homogeneity. 

Higher Multipoles Higher multipoles generate y-distortion of the CMB, as we have seen in Sec. HE. A large 
quadrupole principally indicates large shear, which in turn generates an octopole and higher moments. Using 
the EGS theorem, and the scattering theorem in Sec. IIIB we can see that a small y-distortion indicates strong 
evidence for FLRW geometry under quite general conditions. Scattering of the quadrupole seen by clusters 
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FIG. 14: Without inflation the last- 
scattering surlace consists of 10^ causally 
disconnected regions. We only directly 
observe those on our primary CMB sky 
and we assume the others inside our past 
lightcone share precisely the same history. 
(From [200].) 




Ti = 



Big Bang Singularity Particle Horizon 



induces a polarisation signal in our CMB sky and correlations between clusters can be used to place constraints 
on homogeneity [197]. 

Polarisation The SZ effect also affects the intrinsic polarisation of the CMB [137]. The transverse velocity of the 
cluster relative to the CMB frame can be extracted and used to constrain homogeneity [198] (see also [95]). 

In addition to these features, the CMB can be used with the integrated Sachs- Wolfe effect to constrain the geome- 
try [51, 199], but it is challenging to do so in a model independent manner. 



C. Probing homogeneity in the early universe. 

Homogeneity in the early universe may be constrained in different ways. In particular, the early universe provides 
the primary place we can place constraints on various parameters such as the baryon fraction whose physical origins 
are poorly understood, yet are assumed to be homogeneous. This follows from inflation of course, so probing FLRW 
geometry at early times, and in particular whether widely separated regions share a common thermal history [124], is 
an important consistency check for the standard model. As we have seen, inhomogeneous models require a significant 
reworking of the early universe in order to have a chance of working as dark energy alternatives. Each era provides 
opportunities for confirming homogeneity: 

Last scattering To a first approximation, this happens at fixed temperature on a surface of constant proper time 
in an FLRW model. This temperature depends principally on the baryon fraction and the baryon-photon ratio, 
as well as A^eff) all of which should be constrained to be homogeneous. While we have tight constraints on the 
variability of these parameters across our CMB sky (e.g., the constraints on isocurvature modes are quite tight - 
see e.g., [201]), there are no constraints from within our past lightcone. The temperature of the last scattering 
surface may be probed using the SZ effect, though this must be disentangled from inhomogeneity in the region 
of the scatterer. 

Drag epoch The sound horizon size at the drag epoch also depends on the the (effective) number of relativistic 
degrees of freedom in addition to the baryon-photon ratio and the baryon fraction. This scale is probed inside 
our past lightcone via the BAO. Assuming that any large scale inhomogeneity we are probing has a length scale 
much larger than the horizon size at this time, the proper size of the sound horizon at the drag epoch may be 
written locally as 



dT 



lb 



-1/2 



^ ^ 3 n7fc?7 
4 ijj-fT 



-1/2 



(96) 



where /i/i?o ~ 2998 Mpc The temperature of the drag epoch Td is also a function of fb and rj. This scale is 
imprinted in the galaxy correlation function. Assuming we could measure the Hubble rates i/y and through 
other means, the BAO peak is then a direct measurement of (is(r;e), and therefore constrains fb{x) and r]{x). 
Note that this highlights the problem of using the BAO to infer the Hubble rate at late times: we typically have 
to assume homogeneity at early times to do so in order to find ds- 



BBN Primordial nucleosynthesis is a powerful test of early time homogeneity, because we can measure local light 
element abundances which were synthesised on our past worldline at early times. This can be compared with 
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FIG. 15: B AO and measurements of the CMB 
anisotropies around distant observers tell us 
conditions around the time of decoupling in- 
side our past lightcone, and so allow us to 
constrain homogeneity at that time. To use 
these probes in a model-independent way they 
need to be cross-correlated with other obser- 
vations in different directions such that they 
both measure conditions at the same radius. 
In principle a spatial map of the interior of 
our past lightcone can be made of this era. 
(From [81].) 




spectral distortions of tiie CMB 
or tlie kSZ effect measured 
here is related to the 
temperature of decoupling 



element the element abundance at high redshift and at the CMB - i.e., in the past of worldlines very far from 
us. The ratios of light element abundances depend strongly on the baryon-photon ratio rj which can be probed 
spatially by measuring it in the CMB, and through local abundances of hydrogen, helium and lithium [75, 135]. 
It is intriguing that abundances of ^Li implies a value for the local baryon-photon ratio at odds with the CMB 
observed one by about 50%, and the disagreement is around 5-a. It is easy to interpret this as a spatially varying 
7] [75]. This is one of the few observations potentially at odds with the homogeneity assumption. 

Neutrino decoupling The neutrino background contains information from inside our past lightcone if neutrinos are 
massive. In that case, they do not travel on null geodesies, and therefore contain information about inhomo- 
geneity from the time of neutrino decoupling [47] . 

Primordial power spectrum The amplitude of the power spectrum measured via the CMB compared to low red- 
shift tracers may be different for several reasons. In the standard model, a tilt to the spectrum implies that 
different scales have different power. In a more general model with large-scale inhomogeneity the spectrum in 
one region may be effectively disconnected from the spectrum in another. In the case of a spherically symmetric 
model this is easy to envisage because the spectrum in each shell may have different amplitude, while maintain- 
ing scale-invariance (say) within each shell. Thus, the primary CMB delivers the amplitude of perturbations in 
a comoving sphere ^ 13 Gpc away from us (as well as the tilt, over that sphere), while local measurements of erg 
tell us the amplitude of primordial perturbations smoothed over a comoving region close to our past worldline. 
Tracers at different redshifts probe power in different shells around us, and thus can be used to constrain ho- 
mogeneity at very early times. Inhomogeneity could also be signalled by deviations from scale invariance which 
is different in the radial and angular directions (e.g., a different tilt may be found in the CMB versus LSS), as 
well as from shell to shell. 

The end of inflation/'bang time' Many inhomogeneous models contain a 'bang time' functional degree of freedom 
in them. In a comoving frame such a function has the interpretation of a big bang singularity along each worldline, 
except it happens at different times in the past according to observers on a hypersurface orthogonal to the fluid 
flow. What such a degree of freedom might actually represent in a realistic model remains to be seen, but it 
can be thought of as a 'decaying mode' degree of freedom which can in principle be constrained observationally. 
Recall Eq. (74): Given an observational probe of H^j{z), and a model for local structure evolution (galaxy and 
star formation), this allows us to constrain any bang time inhomogeneity between worldlines A and B [175]. 
Much tighter constraints come from Goodman's tests - it must be 0{t^,) [107, 108]. 



It is the stunning perfection of the CMB that persuades most of us that the FLRW models are the correct basis 
for our cosmological model. Geometrical alternatives which can reproduce such a perfect CMB require fine-tuning 
in space (e.g., LTB models) or time (e.g., Bianchi models [180]), or require ill-defined matter content (e.g., Stephani 
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models [160, 172, 202]). Baycsian arguments suffice to persuade us our working model is most likely correct for 
now. But until we understand A - or whatever dark energy is - with what level of confidence should we ascribe it? 
Rigorously testing the Copernican assumption, and with it the FLRW framework, is therefore a vital aspect of dark 
energy studies. 

We have seen three distinct approaches which help us establish homogeneity: 

Model Building: This has been the most prolific method as it is relatively straightforward to make a model such 
as LTB and compare it to observations. One can then compare an alternative directly to FLRW and decide one 
way or the other. However, this approach has focussed almost exclusively on providing an alternative to dark 
energy rather than testing the Copernican principle per se, so it is not yet clear exactly what conclusions we 
may draw directly. Constraining inhomogeneity when dark energy is present would help test the Copernican 
principle generically (e.g. [151, 203]). 

Consistency Tests: The FLRW models can be formulated as a series of observational null hypotheses which can be 
used to test the Copernican principle (among other things). The idea is to combine observables in such a way 
that they can reveal unambiguously if the FLRW assumption is wrong. Ideally these are independent of dark 
energy or theory of gravity. Where possible, 'observables' such as the distance redshift relation are constructed 
by smoothing data in as model- independent way as possible to remove any residual bias which might exist. 
Robust application of these requires removing as many model dependent assumptions as possible: for example, 
in the case of SNIa this is relatively straightforward; in the case of the BAO this is much harder. 

Inverse Approach: The idea of using data directly to construct the metric on the past lightcone, and then integrating 
into the interior, goes back to a pioneering paper by Ellis et. al. [122] based on a very old thesis by Maartens [125]. 
This is an approach to cosmology which circumvents the Copernican assumption entirely: in theory, with perfect 
observations, we have just enough information from our past lightcone to reconstruct the correct model of the 
universe within it - whatever it is. One of the key requirements is for transverse velocities (or H±ab), so in 
practise this is a very challenging method, but important nevertheless. 

Each of these approaches reveals significant insights into the geometry of the universe, and the restrictions that 
different observables place upon it. However, it is not yet clear from a theoretical point of view when we are able to say 
enough is enough. That is, a principle question of establishing confidence levels remains unanswered. If we constrain, 
say, 'rf{z) < and ^{z) < e^r over some range of redshift, then under what circumstances is that sufficient 
to conclude approximate FLRW geometry? A recent study suggests this might be critical to their efficacy [204]. 
Similarly, if we can constrain a variety of radially inhomogeneous models with dark energy to be homogeneous, is this 
sufficient? 

Some key questions remain: 

Stability: Which of the consistency tests are perturbatively stable? That is, for which of the tests does '= ©(e)' 
imply the geometry is similarly close to FLRW? What does close to FLRW mean? 

Sufficiency: What are the sets of consistency tests which are necessary and sufficient to determine FLRW geometry? 

Confidence: How do we use these to place a measure on model-independent limits of deviations from FLRW, given 
realistic data? 

In the case of exact spherical symmetry and a dust equation of state, say, sufficiency is relatively straightforward ~ 
the LTB model has two free functional degrees of freedom, so if these are chosen to have the same observables as an 
FLRW model with the same matter content, then the models must be the same [7]. This could be phrased in terms of 
the consistency relations as: "if (z) = = J^{z) iff FLRW. Similarly, one can probably say that a model which passes 
Goodmans CMB tests and has "if (z) = is necessarily FLRW - but there is no proof of this as yet. Furthermore, in 
a realistic setting with near isotropy, approximate consistency and unknown dark energy, are these enough? How do 
these tests manage with realistic data? 

It is actually surprising that it is possible to test the Copernican principle at all, and it has only lately been seen 
as a testable assumption. As recently as 2006, George Ellis asserts in relation to the Copernican principle [120]: 
"Establishing a Robertson- Walker geometry for the universe relies on plausible philosophical assumptions. The 
deduction of spatial homogeneity follows not directly from astronomical data, but because we add to the observations 
a philosophical principle that is plausible but untestable." It seems that this is rather too pessimistic. Wc have 
discussed quite general ideas for testing the Copernican principle which can observationally establish homogeneity 
within our horizon which are achievable over the coming years. 
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Covariant formulation of the field equations 

The 1+3 covariant Ehlers-Ellis formalism provides a physically transparent formulation of the field equations in 
fully nonlinear generality (see [205, 206] for reviews). The Ehlers-Ellis formalism is a covariant Lagrangian approach 
to gravitational dynamics, based on a decomposition relative to a chosen 4- velocity field u". The fundamental tensors 
are 

Kb = 9ab + UaUb, Eabc = VabcdU'' , (97) 

where hab projects into the instantaneous rest space of comoving observers, and Sabc is the projection of the spacetime 
alternating tensor rjabcd = — \/— 5^°[a^^b'5^c'^^d]) and so 

Vabcd = "^U^aSbjcd " "^^abicUd] , ^abc^'^''^ = 3\h[/hb^hc]^ . (98) 

The projected symmetric tracefree (PSTF) parts of vectors and rank-2 tensors are 

Via) = K'Vb , Siab) = {hia'h/ - ^ft^'^/labj^cd • (99) 

The skew part of a projected rank-2 tensor is spatially dual to the projected vector, Sa = ^SabcS^'"'\ and then any 
projected rank-2 tensor has the decomposition Sab = \Shab + SabcS"^ + S^ab)- where S = S^dh'^'^- Scalars, projected 
vectors and PSTF tensors are the fundamental objects which represent the gravitational field. 
The covariant derivative Vo defines 1+3 covariant time and spatial derivatives: 

= M^VcJ"-...6, V>,r -...b = h/h'^e ■ ■ ■ hb^VdJ^-...f. (100) 

The projected derivative D^, defines a covariant PSTF divergence, divV — D"V"o , div5o — D^'Sab, and a covariant 
PSTF curl, 

curl Va = Sabc^^V^ , CUrl Sab = Scdia^^Sb/ • (101) 

The relative motion of comoving observers is encoded in the PSTF kinematical quantities: the volume expansion 
rate, 4-acceleration, vorticity and shear, given respectively by 

e = D"Ua, Aa = Ua, UJa = CUllUa, (Taft = D(„M,,) . (102) 

Thus 

VbUa = ^^Kb + ^abci^" + O-ab - ^a^b • (103) 

A key identity (valid in the fully nonlinear case) is 

curl Da V = ea6cD^D'=^/. = -2^/.^^ , (104) 

which shows that curl grad is nonzero in the presence of vorticity (a purely relativistic feature, with no Newtonian 
analogue). A crucial nonlinear commutation relation for scalars is 

C(D(,Vj - = i>Aa - (^Qhab + Oab + Eabe^') ■ (105) 

The PSTF dynamical quantities which describe the sources of the gravitational field are: the (total) energy density 

P = TabU°'u^, isotropic pressure p = ^habT""^, momentum density Qa = —Tia)b'^^^ and anisotropic stress iTab = Tiab)7 
where Tab is the total energy-momentum tensor. The locally free gravitational field, i.e. the part of the spacetime 
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curvature not directly determined locally by dynamical sources, is given by the Weyl tensor Cabcd- This splits into 
the PSTF gravito-electric and gravito-magnetic fields 

Eab = CacbdU''u^ , Hab = ^EacdC^'''beU^ , (106) 

which provide a covariant description of tidal forces and gravitational radiation. 
The Ricci and Bianchi identities, 

'^[a'^b]Uc= RabcdU'^, VCabcd = [a{Rb]c - ^R9b]c} , (107) 

produce the fundamental evolution and constraint equations governing the covariant quantities. Einstein's equations 
are incorporated via the algebraic replacement of the Ricci tensor 

j^ab ^ j^ab _ ^T^Cgab ^ ^gab ^ (^Qg) 

where T"^ is the total energy-momentum tensor. 

The resulting equations, in fully nonlinear form and for a general source of the gravitational field, are: 
Evolution: 

p+{p + p)@ + diyq = -2A''q^-a''^Trab, (109) 
e + + ^{p + 3p) - A - div A = -aabcj"^ + 2WaW" + AaA'' , (110) 

q{a) + ^&qa + {P + P)Aa+'DaP + div'Ka = -(Jabq^ + eabcf^^ - A^T^ab , (HI) 

2 1 ^ 

W(a) + gQ'^a + 2CUrlAa = aafoW , (112) 

2 1 

0'(o5> + gQcTab + -E'ab - ~ ^(o^b) = -<^c{a<^b)'' - '^{a'^b) + At^a^b) , (113) 

E{ab) + - CUrliJ„6 + ]^{p + p)(Jab + \'^{ab) + ^^T^ah + ^D(a96) = -^(a9b> + 2^''£cd(a-f^6)'* 

+ '^'^c{aEb)'' - ^"l^cdiaEb)'^ --j^<y'' {aT^^c- ^'^''^cdia'^b)'^ ^ (114) 
H{ab) + ^Hab + curl Kb - ^CUrlTTab = Sa^i^a^b)'' - UJ^'ScdiaHb/ 

3 1 

-2A''s^a(aEbf - ^l^iaQb) + ^(^^ {a£b)cd<f ■ (US) 

Constraint: 

divw = A'^uja, (116) 
2 . 

divO-a-CUrlWa- -D^e + ga = -2£abc'^ A" , (117) 

curl (Tab + D(aa;6) - Hab = -2A(„a;6) , (118) 

divK + ^divTTa - ^DaP + ^@qa = SabcCT^H'"^ - SH^bUJ^ + ^aabQ^ - ^£a6cW^" , (119) 
diviJa + ^CUrlga - (P + p)0Ja = -eabcT'^dE"' - ]^eabc(j\T^'"^ + ?>EabUJ^ - ^TTabUj'' . (120) 

The energy and momentum conservation equations are the evolution equations (109) and (111). The dynamical 
quantities p,p, qa, Tfai^ in the evolution and constraint equations (109)-(120) are the total quantities, with contributions 
from all dynamically significant particle species. That is, 

j,ab ^ J^Tf ^ pu'^u'' +ph'''' + 2q^''u^'^ +'k"-\ (121) 
Tf = p*ju'iu'i+p*ihf + 2qf"-Ui^ +T^*j''\ (122) 
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where / = r, n, 6, c, A labels the species. The asterisk on the dynamical quantities p}, • • • is intended to emphasize 
that these quantities are measured, not in the u"-frame, but in the /-frame, whose 4-velocity is given by 

W} = 7/ (n« + v^) , v^ua = 0, ^i={l-vj)~'^\ (123) 
The total dynamical quantities in Eqs. (109)-(120), are given by 

p = J2p^^p- . 9" = E«? ' ^""^ = E^r • (124) 
III I 

Assuming that the species are non-interacting, they each separately obey the energy and momentum conservation 
equations (109) and (111): 

P/ + (p/+P/)e + D„g? = -2Aaq] - aabnf , (125) 
9i"^ + ^eg? + (p7+p/)A« + D«p, + D,7rr = -(t V? + e^c^^J - ^bTrf . (126) 

The Ehlers-EUis covariant kinetic theory description starts by splitting the photon 4-momentum as 

p" = Eiu" + e") , e"e„ = 1 , e^Ua = . (127) 

Here E = —UaP^ is the energy and e" = p^'^^ /E is the direction, as measured by a comoving fundamental observer. 
Then the photon distribution function is decomposed into covariant harmonics via the expansion 

f{x,p)= f{x,E,e) = F + F,e" + F,fce"e'' + • • • 

= Y^FA,{x,E)e^^'\ (128) 
e>o 

where Ai = aia2-- -ae and e^' = e""^ ■ ■■ e"''. The multipoles Fa/, are a covariant alternative to the usual expansion in 
spherical harmonics. They are PSTF: 

Fa...b = i^(a...6) ^ Fa-b = -F(a...6), Fa-bU^ = = F^...bch^''. (129) 

The first 3 multipoles determine the radiation energy-momentum tensor. 



Ti^\x) = J p^p''f{x,p)d^p 



= pru'^u'' + \prh"-'' + 2q<fv!''^ + , (130) 

where d^p = EdEdD. is the covariant volume element on the future light cone at event x. It follows that the dynamical 
quantities of the radiation (in the u^-frame) are: 

Pr = An E^FdE, qr = ^ E^F^dE, <'' = —/ E'^F'^^AE. (131) 

We extend these dynamical quantities to all multipole orders by defining the brightness multipoles 

nai...a, = j E^Fa,...a,dE, (132) 



SO that 



U^l-p.,U^^^q^^,U^^='£.f. (133) 



From these the non-linear temperature fluctuations may be defined as, for € > 1, 



Ta,= {^)iIa„ (134) 



where the monopole temperature is defined from pr = aT'^. 
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The coUisionless Boltzmann (or Liouville) equation is 



df 



0, 



dv ^ dx°- dp°- 
where = dx"'/dv, and we neglect polarization. The covariant multipoles of df /dv are given by 



(135) 



1 (dl 
EXdv 
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(2€ + 3) 
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2£ 



{2£ + 3) 



^-1/2 



{Af_2j '^ae-iae) ) 



(136) 



where a prime denotes d/dE. This is a fully nonlinear expression. 

Multiplying Eq. (136) by E^ and integrating over all energies leads to the brightness multipole evolution equations: 



= n 



{t+m 2) ^fcjj^^^ ^ ^ 3) A(„,n^,_,) 



+ 



{2£ + 3) 



{£-m+m+2) , 

(2£ + 3)(2£ + 5) 



(2^ + 3) 



O" {o«n^f_i)6 - (^ + 2)cr(o^a^_^n^^_2) 



(137) 



Once again, this is a fully nonlinear result. The monopole evolution equation is just the energy conservation equation, 
i.e., Eq. (125) with I = r, and the dipolc evolution equation is the momentum conservation equation (126), with 
I = r. The quadrupole evolution is given by Eq. (45). 
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